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^association  along  with  some  entropy  measures  that  provide  diagnostics  for 
|  model  selection  and  testing  for  independence.  Graphical  display  of  the 
estimated  bivariate  dependence  density  and  density-quantile  is  the 
responsibility  of  the  user,  but  a  two  step  procedure  utilizing  the  GCONTOUR 
|  and  G3D  procedures  of  SAS/GRAPH  will  be  suggested  in  a  later  section. 
Different  installations  will  have  different  software  to  provide  three 
dimensional  plots,  hence  the  plotting  option  has  been  excluded  from  BISAM. 

One  of  the  advantages  of  the  modeling  approach  adopted  by  the  BISAM  program 
is  the  ability  to  obtain  quickly  and  efficiently  a  set  of  function  values  for 
a  dense  grid  of  bivariate  points  promoting  the  examination  of  three 
dimensional  representation  of  functions  of  interest. 
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A  USER'S  6UI0E  TO  BISAM:  The  Bivariate 
Data  Modeling  Program 


I.  Introduction 

The  BISAM  (Bivariate  SAMple)  program  is  a  companion  program  to  the 
ONESAM  (Parzen  and  Anderson,  1980)  and  TWOSAM  (Prihoda,  1981)  programs 
developed  for  nonparametric  data  modeling.  The  main  purpose  of  BISAM  is 
to  perform  bivariate  data  analysis  using  Fourier  expansions  and  quantile 
techniques.  The  motivation  and  theory  behind  the  data  modeling  approach 
incorporated  by  BISAM  is  detailed  in  Woodfield  (1982),  with  foundations 
provided  by  Parzen  (1979),  Tartar  and  Kronmal  (1970,  1976),  and  Kimeldorf 
and  Sampson  (1975). 

BISAM  is  a  FORTRAN  program  composed  of  a  main  program  and  37  subprograms. 
Several  of  the  subprograms  come  from  ONESAM  and  the  TIMESBOARD  Time  Series 
Subroutine  Library  (Newton,  1979).  The  main  feature  of  BISAM  is  the  ability 
to  provide  estimates  of  the  bivariate  density-quantile  function  of  a  set 
of  bivariate  data.  A  univariate  analysis  is  also  provided  similar  to  that 
of  ONESAM  but  with  considerably  less  output.  The  univariate  analysis  is 
not  Intended  to  replace  a  ONESAM  analysis  if  such  a  detailed  analysis  is 
required.  The  bivariate  analysis  includes  display  of  various  nonparametric 
measures  of  association  along  with  some  entropy  measures  that  provide 
diagnostics  for  model  selection  and  testing  for  independence.  Graphical 
display  of  the  estimated  bivariate  dependence  density  and  density-quantile 
Is  the  responsibility  of  the  user,  but  a  two  step  procedure  utilizing  the 
GCONTOUR  and  G3D  procedures  of  SAS/GRAPH  will  be  suggested  in  a  later 
section.  Different  installations  will  have  different  software  to  provide 
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three  dimensional  plots,  hence  the  plotting  option  has  been  excluded  from 
BISAM.  One  of  the  advantages  of  the  modeling  approach  adopted  by  the  BISAM 
program  is  the  ability  to  obtain  quickly  and  efficiently  a  set  of  function 
values  for  a  dense  grid  of  bivariate  points  promoting  the  examination  of 
three  dimensional  representation  of  functions  of  interest. 

This  document  will  explain  the  algorithms  employed  in  generating  output 
and  the  steps  necessary  to  perform  an  analysis  of  a  set  of  bivariate  data. 
The  options  available  and  the  user  input  required  to  execute  the  BISAM 
program  will  also  be  described  in  detail.  The  user  should  be  thoroughly 
acquainted  with  the  input  options  before  attempting  to  execute  the  BISAM 
program.  Incorrectly  specified  option  values  may  cause  errors. 

Furthermore,  BISAM  is  a  very  long  and  complex  program  that  will  be  fairly 
expensive  to  use,  and  hence  improper  runs  should  be  avoided. 

2.  Univariate  Analyses 

The  analysis  of  each  of  the  paired  variables  separately  is  essential 
to  fully  understand  the  nature  of  the  bivariate  data  set.  In  particular, 
the  estimation  of  univariate  density-quantile  functions  should  be  performed 
in  an  optimal  manner  to  insure  that  the  estimated  bivariate  density-quantile 
function  is  appropriate.  A  univariate  analysis  will  provide  diagnostics  to 
help  one  understand  the  nature  of  the  univariate  density-quantile  function 
and  thus  guide  one  in  the  formulation  of  the  bivariate  density-quantile. 

BISAM  provides  output  for  each  variable  that  is  a  subset  of  what  may  be 
provided  by  ONESAM.  The  univariate  analysis  is  done  in  two  stages,  however. 
In  the  first  stage,  descriptive  statistics,  an  Informative  Quantile  plot, 
and  a  plot  of  D(u)  are  provided  in  a  goodness-of-fit  exploration  analysis 
to  guide  one  in  determining  the  underlying  marginal  distributions. 


In  the  second  stage,  a  univariate  density-quantile  function  is  estimated 
by  the  autoregressive  method  for  the  null  case  specified  and  is  then  used  to 
form  the  bivariate  density-quantile  function  as  described  in  section  4.  If 
0(u)  in  stage  one  does  not  approximate  a  uniform  distribution  function  for 
the  null  case  specified,  then  one  should  specify  another  null  value,  or  the 
bivariate  density-quantile  function  estimated  will  be  unreliable. 

For  a  complete  description  of  some  of  the  above  concepts,  see  Parzen  (1979) 

3.  Bivariate  Analysis 

A  bivariate  analysis  of  a  set  of  data  includes  computation  of  various 
nonparametric  correlation  coefficients,  model  selection  diagnostics,  and 
estimation  of  the  bivariate  dependence  density  and  density  quantile 
functions  for  a  dense  grid  of  points.  Scatter  plots  are  produced  for 
the  original  data  and  the  rank  transformed  data.  Obtaining  "publication 
quality"  plots  depends  on  the  plotting  hardware  and  software  available  on 
the  system  used. 

Five  correlation  coefficients  are  produced  along  with  various  entropy 
measures  of  association.  We  will  assume  that  the  input  data  set  is 

denoted  by  (X^,  Yj ) . (Xn>  Yn).  Pearson's  product  moment  correlation 

coefficient  is  defined  by 

r  -  l  (X.  -  X)(Y.-Y) 
i*l  1  1 

and  is  labeled  "PEARSON"  in  the  BISAM  output.  Spearman's  rank  correlation 
coefficient  is  given  by 


where  Qi  ■  rank  (X^ )  and  -  rank  (Yj )  and  Is  labeled  "SPEARMAN"  in  the 
BISAM  output.  Three  correlation  coefficients  related  to  the  concepts  of 
concordance  and  discordance  are  also  computed,  two  of  which  perform 
corrections  for  tied  observations.  Kendall rs  Tau  is  defined  by 

tA  =  (NC  '  N0)/[n(n"1)/2:I  . 

where  Nc  is  the  number  of  concordant  pairs  and  Nq  is  the  number  of 
discordant  pairs,  and  is  labeled  "KENDALL  A"  in  the  output.  Kendall's 
TAU  correcting  for  ties  is  given  by 


where  Tx  is  the  number  of  tied  observations  in  the  X  variable  but  not  in 
the  Y  variable,  and  Ty  is  the  number  of  tied  observations  in  the  Y 
variable  but  not  in  the  X  variable.  This  measure  is  labeled  "KENDALL  B" 
in  the  output.  Somer's  d  is  defined  by 

d  =  (Nc  -  Nd)  /  (Nc  +  ND  +  Ty)  , 

and  is  seen  to  be  similar  to  Kendall's  Tau-B.  This  measure  of  association 
is  labeled  "SOMER'S  D"  in  the  output. 

The  ranking  procedure  employed  assigns  average  ranks  for  tied  observations 
Other  methods  for  assigning  ranks  to  tied  observations  are  often  employed  but 
are  not  attempted  by  BISAM.  The  presence  of  a  large  percentage  of  tied 
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observations  will  weaken  the  results  obtained  since  underlying  continuous 
distributions  are  assumed.  One  should  avoid  such  situations  if  possible. 

Two  methods  of  bivariate  density  estimation  are  performed  by  BISAM, 
the  nearest  neighbor  technique  and  the  orthogonal  expansion  technique  of 
Woodfield  (1982).  However,  estimates  obtained  from  the  nearest  neighbor 
technique  are  not  displayed,  but  instead  are  used  as  input  values  for  the 
ortnogonal  expansion  technique.  Since  the  density  estimation  methods 

employed  are  fundamental  to  the  bivariate  analysis  performed,  the  next 
section  will  discuss  these  techniques  with  considerable  attention  paid  to 
the  theory  and  implementation  of  such  methods.  Since  the  entropy  measures 
of  association  and  model  selection  diagnostics  are  so  closely  related  to 
the  density  estimation  technique  employed,  the  discussion  of  these  quantities 
will  also  be  withheld  until  the  next  section. 

4.  Nonparametric  Density  Estimation 

Nonparametri c  density  estimation  is  performed  three  different  ways  by 
BISAM  for  different  purposes.  The  k-nearest  neighbor  technique  is  employed 
first  for  k*6  to  provide  an  estimate  of  the  dependence  density  as  input  to 
an  orthogonal  expansion  density  estimation  technique  and  to  produce  a  raw 
estimate  of  the  entropy  associated  with  the  joint  probability  density 
function  (p.d.f. )  and  the  marginal  p.d.f. 's.  The  orthogonal  expansion 
method  of  Woodfield  (1982)  is  then  employed  to  obtain  smoothed  estimates  of 
the  dependence  density  for  various  degrees  of  smoothing.  Finally,  autoregressive 
estimates  of  the  marginal  p.d.f. 's  are  obtained  and  used  with  the  estimated 
dependence  density  to  produce  an  estimate  of  the  bivariate  density-quantile 
function.  We  now  elaborate  on  the  details  of  these  procedures. 


Let  (X^,  Y-j (Xn,  Yn)  be  a  bivariate  random  sample  form  the  random 
vector  (X,Y)  with  joint  cumulative  d.f.  Fx  y,  marginals  Fx>  Fy,  joint 
p.d.f.  fx  marginal  p.d.f. 's  fx,  fy,  and  quantile  functions  Qx,  Qy.  Defi 

the  dependence  distribution  function  D(u-|,  u2)  by 


D(Uj *  Ug)  -  Fx  y(Qx(u-j),  QY(u2)), 


and  the  dependence  density  function  d(u-|,u2)  by 


d(uru2)  =  3u 


_ D/u  x  -  fX,Y(QX(ul)QY(u2^ 


Define  the  information  between  two  densities  f  and  g  by 
i(f;g)  *  r  Hog  ||*j- }  f(x)  dx 


and  the  entropy  of  a  density  f  by 

H (f )  =  /"  -{log  f (x) }  f(x)  dx. 


The  information  inequality  states  that  for  two  densities  f  and  g. 


I(f;g)  >  0. 


Furthermore,  I(f;g)  *  0  iff  f*g  a.e. 
It  is  easy  to  see  that 


which  justifies  naming  df^^)  the  dependence  density,  since  by  virtue  of 

the  information  inequality  df^^)  is  indirectly  related  to  a  measure  of 

dependence  between  X  and  Y.  This  fact  is  exploited  in  obtaining  diagnostics 

for  model  selection  and  tests  of  independence. 

00 

Let  (0.  (u)}  be  a  complete  orthonormal  system  of  functions  in 
K  k=-» 

o  “ 

L  (0,1).  Then  the  system  {0 .k(u.| .Ug)}  defined  by 

j  ,k=-“> 


®jk(uru2)  =  ^j<ul )<»k(u2)  *  °-ul ,u2-  *  a11  j»kt 


is  a  complete  orthonormal  system  in  the  space  of  bivariate  square  integrable 
functions  on  the  unit  square.  If  log  d(u-|,U2)  is  square  integrable,  then 


log  d(u-|,  u2)  =  2  0jk0j(ul)0k(u2)  "  ’F(!) 

j ,k*-« 


in  the  sense  of  L  norm  where  { are  the  Fourier  coefficients 

JK  j.k— 

defined  by 


ejk  =  fo  fo  ^*WU2*  log  d^ul»u2^  duidu2»  j»k=— »  •••”  • 

The  term  'f(o)  is  included  to  Insure  that  dCu-j.Ug)  integrates  to  one.  For 
the  truncated  m-th  order  model  given  by 
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m 


1o9  "  jjkL„  ®jk>j(ul)0k(u2)  -  ’md)  • 

it  follows  that  log  dm(uj  *u2)-*-  log  d(u^,u2)  as  ik»«,  and  hence  one  calls 
dm(ui,u2)  the  m-th  order  approximation  of  d(u-j,u2).  One  thus  seeks  an 
estimator  of  dm(u]»u2)  based  on  a  random  sample  of  bivariate  data. 

Let 


uli  =  Ri/^n+1)»  u2i  =  (V(n+1) 

where  and  Q..  are  the  ranks  of  X..  and  Y..  as  defined  before.  Thus 

(u-| i •  u2i^*“*»^uin*u2n^  aPPr0X*mates  a  random  sample  distributed  uniformly 

on  the  unit  square  under  an  assumption  of  independence.  Dependence  of  X 

and  Y  suggests  alternate  uniform  distributions  that  may  include  a  variety  of 

bivariate  density  shapes  from  distributions  having  uniform  (0,1)  marginals. 

Next,  form  d(u^,u2)  evaluated  at  these  uniform  sample  points  based  on  the 

k-nearest  neighbor  techniques.  BISAM  uses  k=6.  Form  log  d(u-|,u2)  and  use 

this  as  the  dependent  variable  in  a  least  squares  regression  routine  with 

independent  variables  »  j,k=-m,... ,m.  Essentially,  both 

log  d(U|,u2)  and  (u-| )^k^u2^  rePresent  bivariate  functions  imbedded  into 

univariate  representations  for  a  multiple  regression  analysis.  The  design 

2  2 

matrix  X  is  an  n  x  m  matrix  consisting  of  the  m  orthogonal  combinations 
evaluated  at  the  n  uniform  data  points.  The  routine  performed  by  BISAM 
forms  the  (2hh-1)  x  (2m+l)  correlation  matrix  and  then  uses  a  SWEEP  operator 

♦The  (0,0)  term  is  incorporated  1r  o  the  integration  factor  and  is  not  included, 
i.e. ,  the  no  intercept  m  is  mployed,  hence,  the  correlation  matrix  is  not 
C(2m+1 )2+l]  x  C(2m+1)2  +1J 


on  this  matrix  to  obtain  least  squares  estimates  of  the  Fourier  coefficients. 
The  resulting  estimate  is  given  by 


log  d(u1,u2)  *  ,  j  .  i#*, 


m 

where  {e.  is  the  collection  of  least  squares  estimates  of  the 

j,k=-m 

A  A 

expansion  parameters.  The  integration  factor  4^(0)  is  then  derived  to 
Insure  that  d(uj,u2)  numerically  integrates  to  one.  This  is  accomplished 
for  dtu^tUg)  equally  spaced  in  the  u^  and  u2  directions  within  the  unit 
square. 

Three  models  are  obtained  for  m=l,2,3  which  translates  into  8  variable, 
24  variable,  and  48  variable  regression  models,  i.e. ,  the  terms  with 
subscripts  (1,k),  i,k=-m,...,m  for  m=l,2,3  are  included.  These  three  models 
produce  three  density  estimates  evaluated  at  the  40  x  40  grid  of  (u^ ,u2) 
coordinates  mentioned  above.  Using  raw  Riemann  sums  as  numerical 
Integrals,  one  then  obtains  estimates  of  the  entropy  H(d)  mentioned  earlier. 

A  fourth  estimate  of  the  entropy  of  the  dependence  density  is  provided  by  the 
nearest  neighbor  estimate  using  the  formula 


H(d)  *  /J  /J  lo9  d(uru2)d  Dn(u1,u2) 


where  0n(u1,u2)  is  the  empirical  c.d.f.  of  (u^,  u^),  i«l,...,n. 


A  A 

As  a  rule-of-thumb  for  model  selection,  H(d)  is  compared  with  H(dg), 

A  A 

H ( d24 ) >  and  H(d4g).  An  information  criterion  modeled  after  Akaike's 
Information  Criterion  (AIC)  familiar  to  time  series  analysts  is  used  to 
select  the  "best"  model.  One  forms 

AIC  (k)  =  H(dk)  -  H(d)  -  2k/n 

and  then  selects  model  k  for  the  k  value  at  which  AIC(k)  achieves  its  positive 
minimum.  An  alternative  approach  chooses  k  for  which  |AIC(k)j  is  minimum. 
Model  selection  criterion  are  still  under  investigation  with  promising  results 
having  been  achieved  for  the  criterion  employed  by  BISAM. 

Graphical  displays  of  dk  suffer  from  the  unruly  behavior  of  dk  near 
the  boundary  of  the  unit  square,  and  hence  one  may  prefer  to  display 
values  of  the  bivariate  density-quantile  function  given  by 

f Q ( u1 , u2 )  =  dk(ultu2)  fQx(Uj)  fQY(u2) 

where  fQx  and  fQy  are  the  autoregressive  estimates  of  the  univariate 
density-quantile  functions.  Figure  1  contains  a  three  dimensional  plot 
of  the  bivariate  density-quantile  function  of  order  8  for  a  set  of 
simulated  bivariate  normal  random  variables  with  correlation  coefficient 
p  *  0.9.  Figure  2  contains  the  theoretical  bivariate  density-quantile 
function  for  the  simulated  data.  The  plots  were  produced  using  SAS/GRAPH. 
Appendix  B  outlines  a  two  step  FORTRAN-SAS  procedure  that  allows  one  to 
produce  three  dimensional  or  contour  plots  of  bivariate  functions  of 


Interest. 


5.  Using  BISAM  to  Analyze  Time  Series  Data 

Although  primarily  intended  to  analyze  bivariate  data  in  which  the  data 
values  are  uncorrelated.  BISAM  can  also  be  used  to  analyze  time  series  data. 
For  a  univariate  data  set  one  may  use  option  IUNIV  to  do  an  analysis  of  a 
bivariate  data  set  consisting  of  the  univariate  data  set  paired  with  lagged 
values  of  the  data  set. 

Consider  a  Guassian  time  series  X(t)  with  autocorrelation  function  p(v) 
given  by 


p(v)  =  Corr[X(t),  X(t+v)3,  v=0,  +1 . 

Define  the  series  Y ( t)  to  be  a  lagged  version  of  X(t),  i.e.,  Y(t)  =  X(t+v). 
Then 


p„  M  =  Corr[X,Y]  =  Corr[X(t),  X(t+v)]  =  P(v). 

For  a  bivariate  normal  random  vector  (X.Y)  with  correlation  coefficient 
Pv  it  can  be  shown  that 

'<f*.y  •  Vy>  *  -05  l09  (,-‘>S.y)  • 

Consequently,  the  information  between  X  and  lagged  values  of  X  is 

I(v)  =  -0.5  log  (1-p2(v)) 


In  the  non-normal  case,  this  equation  does  not  hold  so  that  a  nonparametric 


estimator  may  be  desired  to  indicate  the  relationship  between  X  and  lagged 
values  of  X. 

When  one  specifies  1UNIV  =  1,  values  for  KLAG1  and  KLAG2  must  also 
be  specified.  Then  a  bisam  analysis  is  performed  on  (X(t),  X(t+v))  for 
integer  values  v  satisfying  KLAG1  <  v  <_  KLAG2. 

6.  Input  Options 

The  following  options  are  input  on  the  first  data  card  in  1115  format. 
They  are  input  in  the  order  listed,  and  if  NT APE  =  5,  the  data  set  follows 
this  card  in  the  indicated  format  listed  at  the  end  of  this  section. 

NTAPE  -  number  referring  to  DD  statement  describing  the  input  data 
set. 

IDQX  -  null  distribution  for  autoregressive  smoothing  for  X  input 
variable. 

IDQY  -  same  as  IDQX  except  for  Y  input  variable. 

MORD  -  maximum  autoregressive  order  to  be  used  for  univariate 

autoregressive  density  estimation  (<_  6). 

IPLT1  -  0  for  no  scatter  plots. 

1  for  scatter  plot  of  data. 

2  for  scatter  plot  of  rank  transformed  data. 

3  for  both  scatter  plots. 

IPLT2  -  0  for  no  univariate  density  plots. 

1  for  best  order  AR  univariate  density  plots. 

0  for  no  univariate  descriptive  statistic. 

1  for  descriptive  statistics  for  both  X  and  Y. 
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KDEL  - 


IOUTD  - 


I  REG  - 


IUNIV  - 


maximum  number  of  extreme  points  to  exclude  from  bivariate 
analysis.  An  extreme  point  is  located  based  on  the  distance 
of  the  X  or  Y  coordinate  from  its  median.  If  KDEL=2  and 
both  the  extreme  X  and  Y  are  paired  together,  then  only  that 
point  will  be  omitted.  If  KDELs3,  two  points  with  an 
extreme  X  value  and  one  point  with  an  extreme  Y  value  will 
be  omitted,  i.e.,  the  X-direction  received  precedence  over 
the  Y-di recti  on  for  odd  values  of  KDEL.  This  method  of 
deleting  outliers  works  well  for  nearly  linear  relationships 
between  X  and  Y.  Its  use  is  questionable  for  other  cases. 

0  if  function  values  are  not  to  be  saved. 

1  if  function  values  for  the  three  fitted  models  are  to  be 
written  to  tapes  1,  2,  and  3  respectively. 

If  10UTD=1 ,  the  JCL  must  contain  three  DD  cards  for  FT01F001, 
FT02F001,  and  FT03F001  defining  permanent  dish  data  files 
where  the  function  values  are  to  be  written.  Typically,  10 
tracks  of  storage  must  be  alloted  for  each  file. 

0  if  no  quantile  regression  performed. 

1  if  nonparametric  estimates  of  rQ(u)  =  E[Y|X  =  Q(u)]  and 
r(x)  *  E[Y|X=x]  are  desired. 

0  if  input  data  set  is  bivariate  data. 

1  if  input  data  set  is  a  univariate  time  series  as  described 
in  section  5. 


If  IUNIV=1,  a  second  data  card  is  required  containing  the  values  KLAG1  and 
KLAG2  in  215  format  as  described  in  section  5. 


If  IDQX  and/or  IDQY  specify  the  Pareto  to  Weibull  distributions,  a  separate 
data  card  must  be  inserted  containing  estimates  of  the  parametrs  in  those 
distributions.  For  the  Weibull,  the  location  estimate  is  given  first,  with 
the  input  format  specified  as  2F10.0. 

The  input  data  set  consists  of  two  univariate  data  sets  containing  the 
same  number  of  values.  It  is  always  assumed  that  the  X  data  set  is  read  in 
first.  The  values  ( X , Y)  are  then  paired  by  the  order  they  are  read  in.  A 
univariate  data  set  consists  of  a  title  card  read  in  20A4  format,  a  description 
card  containing  the  number  of  data  points  and  the  format  of  the  input  data 
read  in  (15,  4X,  5A4)  format,  and  the  data  cards  containing  the  data  coded 
in  the  format  indicated  on  card  3.  If  the  number  of  data  points  indicated 
for  X  and  Y  do  not  match,  theprogram  will  terminate. 

Example  of  univariate  data  set: 

YEARLY  SNOWFALL  IN  BUFFALO,  1910-1972 
63  (6F10.2) 

71.2  69.5  47.8  58.4  29.9  42.5 

etc. 

The  code  for  BISAM  and  sample  output  are  contained  in  the  Appendices. 

7.  Sample  Output  form  BISAM 

Following  is  a  listing  of  the  output  from  BISAM  for  a  typical  run  with 
the  input  options  clearly  labeled.  The  JCL  for  executing  BISAM  at  Texas 
A&M  University  is  given  in  Appendix  A. 
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APPENDIX  A  -  JCL  for  Executing  BISAM 


//  Job  Card 

//*  JES3  Control  Cards 

//PROCLIB  DD  DSNsUSR .  R579 . TW . P ROC L I B ,  D I S P=SHR 
//  EXEC  BISAM 
//SYSIN  DD  * 

Parameter  Input  Card(s) 

Data  if  NTAPE=5 

If  NTAPEj^,  the  value  of  NTAPE  Is  coded  as  a  number  nn  between  8  and  99 
except  for  11.  (The  value  11  is  currently  used  for  a  scratch  file.)  A  DD 
card  is  then  inserted  defining  the  input  data  set.  For  example, 

//FT13F001  DD  DSN=USR. R579.TW. DATA.DISP-SHR 

described  the  dish  file  containing  the  input  data  with  NTAPE=13.  If 
IOUTD-1,  three  such  cards  will  be  required  as  described  in  section  6. 


£  * 


,  ^  *-7  rz  ^  «-,  •  .  . 


APPENDIX  B  -  OBTAINING  THREE  DIMENSIONAL  PLOTS 


//  JOB  CARD 

//*  JES3  CONTROL  CARDS 

//PROCLIB  DD  DSN=USR.R579.TW.PROCLJB,DISP=SHR 

//  EXEC  BISAM 

//SYS IN  DD  * 


***  OPTION  CARD  FOLLOWED  BY  DATA  IF  NT APE  =  5  *## 
//FTOlFOOl  DD  DSN=WYL.XX.YYY. STOR1 >DISF'=OLD 
//FT02F001  DD  DSN=WYL . XX . YYY . ST0R2 , DI SP=OLD 
//FT03F001  DD  DSN=WYL  .  XX  .  YYY  .  STOR3  *  D I  SP=*OLD 
//FT13F001  DD  DSN=WYL . XX  *  YYY . DAT1 t DI SP-OLD 
/* 

//STEP2  EXEC  SAS »  COND®  <  0 » LT ) 

//FILE1  DD  DSN=WYL,XX.YYY.STORl»DISP=OLD 
//FILE2  DD  DSN=WYL.XX.YYY.  ST0R2  , DI SP-OLD 
//FILE3  DD  DSN=WYL .XX  » YYY • ST0R3 r DI SP-OLD 
//SYSIN  DD  * 

DATA  ONE »  INFILE  FILE1J  INPUT  U1  U2  DOHT  DHAT i 
TITLE  . F=COMPLEX  ,H=1  FIRST  ORDER  BIVARIATE  DENSITY r 
PROC  63D  DAT A=ONE  GOUT  =  A»  PLOT  U1*U2==DQHTJ  . 

PROC  G3D  DAT A=ONE  GOUT=B»  PLOT  U1*U2=DHAT» 

DATA  TWO  f  INFILE  FILE2J  INPUT  U1  U2  DQHT  DHAT i 

TITLE  .  F=COMPLEX  .H  =  l  SECOND  ORDER  BIVARIATE  DENSITY  » 

PROC  G3D  DATA=TWO  GOUT=C5  PLOT  Ul*U2=DQHTi 

PROC  G3D  DATA=TWO  GOUT=D>  PLOT  U1*U2=DHATJ 

DATA  THREE  5  INFILE  FILE3J  INPUT  U1  U2  DQHT  DHAT » 

TITLE  ♦ F “COMPLEX  ,H=1  THIRD  ORDER  BIVARIATE  DENSITY » 
PROC  G3D  DATA=THREE  GOUT  =  E i  PLOT  U1*U2=DGHTJ 
PROC  G3D  DATA=THREE  GOUT=F>  PLOT  U1*U2=DHAT? 

DATA  COMBINE;  SET  A  8  C  D  E  Fi 
PROC  GREPLAY  DAT A=COMBINE » 

/* 

NOTES?  1.  USE  CURRENT  SAS/GRAPH  JCL  IN  LINE  16* 


3.  XX. YYY  IS  USER'S  ACCOUNT. 


5.  THE  SAS  PROCEDURE  GCONTOUR  MAY  BE  SUBSTITUTED  FOR  G3D  TO 
PRODUCE  CONTOUR  PLOTS  INSTEAD  OF  3D  PLOTS. 

6.  AS  INDICATED »  PLOTS  OF  DHAT  MAY  NOT  BE  INFORMATIVE  t  SO 
DHAT  PLOTS  INDICATED  ABOVE  ARE  USUALLY  OMITTED. 

7.  IN  LINES  11-13,  THE  WYLBUR  FILES  ARE  DUMMY  FILES  SAVED  IN  CARD 
IMAGE  THAT  WILL  BE  WRITTEN  TO.  IN  LINE  Ml,  NTAPE-13  HAS  BEEN 
SPECIFIED  AND  0AT1  CONTAINS  THE  BIVARIATE  DATA  SET  TO  BE 
ANALYZED. 
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C  • ISAM  ■  1  • AM 

C  HIM 

{••••••••••It . . . . . .  *!•  AM 

c  a i a am 

c  baivir  faobram  m  bivaaiati  bata  anal vs  i  b  bisam 

c  ■  I  BAM 

C  INFUTt  NT  AB|  *  TAM  WHIM  BATA  BIT  MB  I  BIS  BIBAM 

C  |S.V|  •  StVARlATI  BATA  (  I  NO  1  V  IBUAVL  V  .  I  UNIT)  BIBAM 

C  MBA  0  •  MAI  I  MUM  AUTBAICRIBBIVI  0 AOI A  TO  Bl  MB IB  f BA  BIBAM 

C  UN  I  VAR I AT|  AA  OINBITV  IBTIMATIBN  (<<4J  BIBAM 

C  IBOI. IBOV  •  NULL  D I BTA 1  BUT | BNB  9  BA  AUTBAI6AIBB 1 VI  SMOOTH  I  NO  BIBAM 

C  INLTt  -  B->  NO  BCATTIA  BLOTS  BIBAM 

C  1*->  BCATTIA  BLOT  BB  OATA  BIBAM 

C  f**>  BCATTIA  BLOT  OB  BANK  TRANS BO AMID  BATA  BIBAM 

C  *-->  BOTH  BCATTIA  BLOTS  0  ISAM 

C  IBLT*  -  #••>  NO  AUTBAISRIBS I V|  BINBITT  BLOTS  BIBAM 

C  1 • *  >  BIST  BABB A  AN  BINBITT  BLOTS  SIBAM 

C  1 BB T  *  A-->  NO  UN  1  VAR 1 AT I  OIBCA1BT1VI  STATISTICS  BIBAM 

C  (•->  UN  I  VAN  I  ATI  RISC A  1  NT  I  VI  STATISTICS  BOR  I  AND  T  BIBAM 

C  NOIL  *  MAI  I  MUM  NUMB I A  OB  I ■ T  A  C  M|  BOINTS  TO  BKCLUDI  BROM  BIBAM 

C  B 1  VAR  I  AT |  ANALVSIB  BIBAM 

C  I  OUT  B  -  t  |B  TNI  S  MOBIL  B  ISTIMATCD  A Al  TO  NAVI  VALUES  BIBAM 

C  MR  t  T  TIN  TO  TABIS  1.1,  AND  1;  •  0 . *  BIBAM 

C  FTOIOOI.  I T  C  .  BO  UCL  CAROS  MU8T  01  INCLUDIO  IB  SIBAM 

C  ISUT0>1.  BIBAM 

C  IRIS  •  1  IB  OUANTILI  ABCACBSIBN  BIRBBAMIO.  O  DA  BIBAM 

C  I  UN  I V  •  t  IB  UN  I  VAR  I  ATI  OATA  SIT  TO  Bl  RIAD  IN  AND  A  BIBAM 

C  B I  V A A  I  AT |  BATA  BIT  CRIATIO  AS  B A  I S INAL  OATA  0  I  BAM 

C  AND  LA66KB  BATA.  •  O.W.  BIBAM 

C  KLACt.Ai.A42  *  R ANSI  BB  LABS  INNUT  IN  I  UN  I  V • I  0  ISAM 

C  BIBAM 

C  B URN ABB RAMS  CAL  L IB  :  ARIST , AUTSON . AUTBRS . CLOLT 1 . CMB INB , CBLOT .  BIBAM 

C  CBTINT . CBORIO . CSwBBF .BATAIN. BISTAT .NCBOBA.  BIBAM 

C  BOA  I IR . BOBNC . BTIAN .  ICOBIA , NINOAL  . ABO . MAI .  BIBAM 

C  MONAIB.MIN.MINMA1 . OAB1. BAA* . BIAABN .BLOTIV .  BIBAM 

C  BBLIT .BBINO.BBLBT .ORIS, OTOBB.OUBNT.OUICN.  BIBAM 

C  RANK.  SNRMN.  TRIM.  NTSBACI  BIBAM 

C  BIBAM 

C  NOTI :  SCRATCH  TAB!  NUMRIR  IS  SIT  AT  11  IN  BCOOIA.  A  00  CARS  BIBAM 

C  BOA  BT11BOOI  MUST  01  INCLUDIO  IN  TN|  JCL  SIT  UB .  BIBAM 

C  BIBAM 

C  UULV  IBOS  BISAM 

C  NA06AAMMIR:  TIRRV  J.  WOOBBIILB  BIBAM 

C  BIBAM 

. . BIBAM 

C  SIBAM 

COMMON  /SATAN/  A  I  BOB » . V | BOB I . BANNS | BOO ) . RANK V  I  BOO )  BIBAM 

COMMON  /BARM/  BITAB. BITAW  0  I  BAM 

0 1 MINS  I  ON  L|S|. LAST | BO )  , L ABA  I BB )  . ▼ I  BOO .2),M0t4|,AIC(4)  BISAM 

0  I M|NS I  BN  CHARIS) . NNAMI | VB) . VNAMI | B0|  . CABTI IB|  . CRNR(B)  BIBAM 

DIMENSION  Ml  'OOO)  . LABOI 2)  . LABOR | 20) . LABOT I  BO)  BIBAM 

IBUIVALBNCI  ITII.II.IID)  BISAM 

•OUIVALINCI  ITfl. S). VIII)  BIBAM 

BATA  NOUT .NIN/S.S/  BIBAM 

DATA  LABB/4M  A  ♦.AN  LAS/  BIBAM 

OATA  CHAR/ IN* . 1N« , 1HA. 1N». IN . /  BIBAM 

BATA  AN AMI / IB* IN  .1HR.B*1M  /  BIBAM 

BATA  VNAMI / 1 #• 1 N  . INV.0*1H  /  BIBAM 

BATA  LABDA/AN0RI6. ANlNAL . AN  BAT . ANA  •  .ANA  .  1 B*  AN  /  BIBAM 

BATA  LAB0V/AN0RI6. ANIMAL .AN  BAT . ANA  •  .ANT  .  1 B*  AN  /  BIBAM 

BATA  CABT/ANBCAT. ANTBR  , ANBLST , AM  OB  .ANA  VB , AN .  V  .  1 A • AN  /  BIBAM 

OATA  CRNK/AN •  AA.AMNK  T  . A NR ANB . ANN 0 AM . ANA T I  0 . A NN  /  BIBAM 

C  SIBAM 

LABTaB  BIBAM 

MR  I  Til NOUT . 1  )  BIBAM 

1  FORMATI INI |  BISAM 

RIABIMIN ,  IB |  NT AOI ,  I  BOA .  I  00 V . MBRO .lOLTI.IOlT*.  I  BBT , NOIL .  IOUTB,  BIBAM 

♦  I RIC , I UN  I  V  BIBAM 

10  9 BRMAT |  1 B I S  |  BIBAM 

IB  I  IUNI V .10 . 1  )  Rl AO  ININ. lO)  SLAB  1  . NL A6*  BIBAM 

t  B |  IBBA .10. t  )  Rl AB ( N I N .  12)  BITAO A  BISAM 

|0t  IBBI .10.  IB)  RIAD(NIN.ia)  BITAOI , BIT AMI  BIBAM 

IBI IDOV.IO.B)  RIABININ.IB)  BITAO V  BISAM 

|0 | IBOV . 10 . 10>  RIABININ.IB)  BIT AO V . Bl T AMV  SIBAM 

IB  BORMATHO  IO  O»  SIBAM 

MR  ITHNBUT.IOl  BISAM 

20  FORMAT | //, lOA . I0| AN* ■■•)./. 10A 7BA .*•*,/, IO*. '•  BIBAM  *,  BISAM 

♦  ' -  BIVARIATI  BATA  ANALVSIB  USING  B  OUR  I IR  I  A BANS  I BNB '.  I  .  BISAM 

♦/.ION.**  ANB  OUANTILI  TICNN lOUIB’.AAA, *  •  * ,/ , IOA . ' • * ,  BISAM 

♦  VBA.  *•*,/.  ISA.  *•  0RB6RAMMIR:  TIRRV  J.  NOB OB  1 1 L  0 '  .TOO.  *  • ' .  BISAM 

♦/. 108.  *• * . TBS.  ••*,/, 108, 20  I  AN* • • • ) I  BISAM 

|0f lUNIV.IO.O)  SB  TO  20  BISAM 

CALL  BATA I N I  NT AO I . W. NM, L . LABS )  BISAM 

BB  200  SLAB  *  KL AS  1 . RL  AGS  BISAM 

N*NM«NLA8  BISAM 

BB  22  1*1. N  BISAM 

Rll)aM(l)  BISAM 

V(I)*M|| *RL  AS )  BISAM 

22  CBNT | NUS  BISAM 

BO  a*  1*1.20  BISAM 

2A  L  AS  V |  I  ) *  L ABA (  I  )  BISAM 

CALL  I CBBIA ( SLAB . AH |  I A )  .LABV|1T)|  BISAM 

LABVI  IB ) *LABB|  1  I  BIBAM 

LABV(IS)*  LABS  12)  SIBAM 

SB  TO  «•  BISAM 

20  CBNT I NUB  BISAM 

CALL  BATAININTAOI. I.NR. L.LABS)  BISAM 

CALL  BATAININTAOI. V.NV.L. LABV)  BISAM 

N«NS  BISAM 

1 0 ( NS  BO.NV)  SB  TO  40  BISAM 

WRITS (ROUT . 20)  LABI. LABV  BIBAM 

20  FORMAT |  IN  ,  1BA.2BA4/.  1OS.20A4.//. lOS.  'SAMOLC  S I  IBB  NOT  IOUAL .  * ,  BISAM 

«'  BIVARIATI  ANALYSIS  I NAORBOR I ATI .  ISICUTtON  TIRMt RATIO .  *  I  BISAM 

B  TOO  BISAM 

40  If  I  LAST .61 . 1 !  SB  TO  St  BISAM 

WRITS | NOUT . BO)  LABI . L AB V . N . N T AOI .  I  BBS. I OS V. MONO, I0LT1.I0LT2.  BISAM 

♦  I BBT .NOIL . IBUTO . IRIS . IUNI V  BISAM 

SO  FORMAT!///. 10A.20A4,/. 108.2004,//. 101. ‘SAMOLI  Bill  •*,IB,  BISAM 

♦  //. 101.  'RATIONS  FOR  THIS  ANAL  VS  1 S  :  ' .  /  .  T»0  ,  '  NT  AM  a  *.12.  BIBAM 

« TOO ,  *  I  BOR  •  ’. 12. TOO.  * IBOV  a  *.12./, T20 ,  *  MBRB  •  *.11.  BIBAM 

♦  TOO,  *  I OL T  1  a  * . I  2 . TOO, * IOi T2  •  * .  I  * . / . T SO .  *  I BST  •  *.12.  BISAM 

♦  TOO,  *ROf L  a  * . 1 2 . T SO .  *  I OUT B  •  * . 1 2 . / . T f O ,  *  I  RIB  •  '.12. TOO.  BIBAM 

•* IUNIV  •  * , I2|  BISAM 

If (I I0LT1 .10. 1  I . BB .11 OLT I .10.2) )  BISAM 

•CALL  OOLBTIS. V. SOB. N. 1. CHAR, CAOT.INAMt, VNAMI. 0)  BISAM 

Bl  WRITIINBUT. 1 |  BIBAM 

C  BISAM 

C  BROS A  BIVARIATI  BATA  W  I  VALUIS  BIBAM 

C  BIBAM 

CALL  BRB2IT.N.SOO)  BIBAM 

If* | IDST .10.0)  .OR. HAST. SI. 1  I)  SB  TO  Bl  BISAM 

BIT AOa  BIT AOS  BIBAM 


SAM 


CALL  OUCNT I M . I . I AIM , LABOR ,  I  001 . M*4. X>S. BMIO . 178 .ISAM. Sil | 


HIM  HT*  BIT  OP  AT  MBIT  Mil  llTMMl  POINTS 
•  >  CALL  TM I M| R , T ,  IMBI,  TMBB  .  ROB  l  ,  M  ,  NIWM  } 


CAL  l  UHKU.V.IAMIt 
CALI  A AMR  I  T  ,  M  , MAMMY  I 

COMPUTB  CBBRBLATIBM  CIlfMCIINtl 

CALL  SPBMMIN. AM0.8UMNI 

CALL  RINOALIM.tAUA.TAUl.BOWSA.R 

CALL  PIAASNlN.A) 


If 
•  • 
A© 


OBTAIN  BBTIMATBB  99  B 1  TAM  I  ATI  MB  Ml  MBBMCB  BfMBMT  AMP  0BNB1TT 
•  UAMTILC  9 UMC T I  BN 

CALL  CMP  I NP I M. MONO.  IBpl.  I  00  V .  I  ALTS.  I  BUT  B ,  I  MB 6 . 1 ABB . L APT . MB . A I C \ 
WAITIIMOUT,  I  > 

IPI(NINOBOO)  AMB. IMBIP.SO.©) I  SB  TB  BB 
WN|TI|NOUT .SB)  NINO, NOIP 

P  BMMAT I  /  /  . '•» , * T IBB  !M  I  T IBB  >M  T  •*.»4.//> 

WM  ITBIMBUT  ,  BB  )  L  All ,  L  AB  V  ,  M 

f BMMAT I//. tOM.SOAA,/. IOI.SOA4.//, IBS. 'SAMPLt  Bill  •'.!§) 

BBT AIM  I  MM  BMMAT I  ON  MEASURES  MOM  MOMMA L  CABC 

N  t ■ *0 . B* AL 06  I  I . -A-Al 
AMB I • *O  I* ALOB  I  I  .  • MMO*MNO I 
T AUA1 » • O . l> AL06 t  \  .  • T AM A* T  AUA I 
T  AUO 1 •  *O  t*ALOCM.  - T  AUO • T  AUB  I 
SBMBMI •  -0.0 •ALOB (  t  .  • B  AMI M • t  BMf  A ) 


:  ii 


MAI TB I  MOOT A . Ml . MHO . AMB I  . T AUA, TAMA I  . TAUB , T AUB 1  . BOMB A , BOMBA  I 
9 BMMAT | //.TIB.  *  C SAAB L AT  I BM  COS? ICIBNT' .TAB,  *  TALUS  *  , TBB , 

«* IMPBRMAT1NM  (NORMAL  CABI I • , / ,T IB. BB| TN* ) , TAB, 1B| IN* ) , TBB . 

♦  SSI  IN- |  . / . TIB.  *  PI ARSON ' , TAB . P IB . 0 , TBB. P 10 • 0 , / . TIB.  *  BAB AMMAN* . 

♦  TBB. 9  IB. B, TBB. 9  IB. 8./. TBB.  ‘RANDALL  A*  , TAN, Pi©. 4. T  AB , MB . I, / , 

♦  TSB.  ‘RBMBALL  B *  . T AB . 9  IB . A . T AB . 9% B . B . / , TBB , B NBOMflA * B  B.TAB, 


T8T, *A!C* . / . T I B , B I  IN* I, 
TIL  BA '.TSB. Alt. S. TAB. MO. B./( 

B • MAT  14‘ , TJS.P1I  .S. 


•P1B. S. TSB. 

P 

IS 

41 

BO  SS  1 ■  1  . 

4 

OS 

MB I  1  1 • *  MO  1 

1 

1 

WRITBIMOUT 

AS 

J  MB  I  A 

1  .  A 

ICI  •  > 

.MO 

MI.AI 

ICI  1  ) 

♦A1C1S) 

SB 

PBAMATI //. 

T 

<B 

.  *  MOBB 

L  ‘  . 

TSB  .  * 

INP 

OAMAT! 

IBM*  . 

•  TSB .  1 1  I  IM- 

1 

. TAD. IBI 

IM* 

l./.T 

ts . 

'B-TILBA*  . 

♦TIB. 'B-MAT 

B* 

.  T SS  .  P 

1  1  . 

S  .  TSO 

.  P  1 

B.S. /. 

,TIO. 

♦TSB. P IB. S. 

/ 

.  T 

10.  *  B * 

MAT 

48  *  . 

T3S 

.Pit  S.TSO 

BB  AB  !•«. 

N 

If  1  I ■RANRI 

1 

1  ) 

/PL  BAT 

IN* 

1  ) 

T|l  ) • AANR  T 

1 

1  I 

/PL  OAT 

IN* 

»  1 

BB 

OB 


CBNTIMUI 

BB  OB  1 • I . A 

CAP T |  I • A  I • CAMS  I  I | 

CART  1  NUB 

IP I  I IPLT1 .10.11 .BA.  I  I  PIT  I .80.31  I 


Ml 


LAST* 1 
lAITHA 

CANT | HUB 

BTBP 

BMB 

PUMCT1BM 


PUMCTIBM  TO  COMPUTB  AMT BMC BMB SB  I TB  BBTIMATBB  BVALUATIO  AT  A 
MBTMBB:  A  A  TWO  ■  OPTRMM  /  ABSM  ♦  ▼)••> 

WNSAt  T  •  APT  CAB lJ)»IXP|  l*J»l*PJ»R)  BUMMBB  B VBA  J  ■  1.  L 

INPUT; 

B:  8CALAA  AT  WHlCM  AUT0AISAB8S I 71  BBT IMATB  IS  E7ALUATI0. 

L:  BABBA.  MUST  BB  LBBS  THAR  II.  SBB  MBTMBB. 

BPTANM:  SKI  MSTWOO . 

BP T C 01 :  AUTpNBSABBBIVfl  COBPPICI BMT8  BP  BABBA  L.  SBB  MBTMBB. 
BPTCB8  IB  A  COMPLEX  10-7BCT0A. 

OUTPUT;  PUNCTIBN  MB  TURNS  TALUS  BP  AUTBABCABSS I TB  BSTIMATOA  BVALUA 
AT  A. 


COMPLIR  bptcbbili 

CAMP LSI  6 
P I  •  A . • AT AM | I .01 
6'CMPiB | I . ,B. | 

BB  I  JH.l 
PJ»4 

6*S*0PTC0II Jl •CBBPtCMPLR (O. , l>2 . >PI*PJ| | 
I  CBNTIMUI 

BBIBT • BP TR MM/ RIAL  I 8a COM JC I C )  I 
IBB 


BM  BASIS  BN  TNI  AUTORIC MISS  I TB  MBTMBB  OP  PAASIM |  I07B I . 

W  -  MAW  BATA 
M  -  SAMPLI  Bill 

IOON  *  f  MB  I C ATOM  PSA  NULL  BIST.  BP  W 

IBBAT  ■  B  IP  W  AMB  AANRW  80MTB8 ,  1  OTMIRWtSI. 

MO AO  •  MAXIMUM  ALLBWABLB  OROIA  (081 
IPLTS  -  O- - >  MB  PLOTS 

1  • ->  PLOT  BP  AM  BINS  I TT-OUART I  LI  PUMCTIBM 
WLAB  •  TAN  I ASL I  NAMB  PON  W  IM  AS  P  BMMAT 

:  MTW  -  BABIN  BP  AUTBABCABSS I T|  BINS  ITT  BSTIMATOA 

A L  PM  *  CBIPP I C I INTB  PON  AUTO All NIB 8  I TB  AB PAB BINT AT  I  OB 
AT A AW  •  AIStBUAL  TAAIANCI  POA  BIST  BABIN 

BIBB  -  I RTICAAT I  MB  PACTBA  I B I 6MA  *  T 1 L  BA  PBA  NULL  MOBIL) 


U| 


COMMON  /PANM/  BITAP, BITAW 

BINIMBI  BA  WIN) .NOANIB) .UISBB) .BNIBBO) .BLIBBB) .PBIBBBl . 
PfOSBIBBB  )  ,  CW88 18001 .  ILBC  IB  I .  CAT  I B  )  .  WR  1 1  BOB  I 
BIMINI I BM  CAPTISO) 

COMPLII  All) ,PNI fSI.ALPHlB) . ALPMAf «B) .NIB VAN 
BATA  C APT/AMUN IT, OMAN  I A , ANTI  B . AMINS I , BMTT -B . ANUANT . AN I L I 
POMPON  , AMNANB.ANBM  V , OMAN  I  A . 4 MB LI  . B*  AM  / 

BATA  SPCPAC/B , I / 


ISAM 

ISAM 

ISAM 

ISAM 

IBAM 


ISAM 

IBAM 

IBAM 

ISAM 


IBAM 

IBAM 

IBAM 

IBAM 


IBAM 

IBAM 


IBAM 

IBAM 

IBAM 

ISAM 

IBAM 

IBAM 

ISAM 


ISAM 

ISAM 

ISAM 

ISAM 

ISAM 


ISAM 

IBAM 

ISAM 


ISAM 

ISAM 

ISAM 

ISAM 

ISAM 

ISAM 

ISAM 

ISAM 

ISAM 

ISAM 


ISAM 

ISAM 

ISAM 

IBAM 


ISAM 

ISAM 

ISAM 

IBAM 

IBAM 

ISAM 

IBAM 

IBAM 

IBAM 

IBAM 

IBAM 

IBAM 

IBAM 

IBAM 

ISAM 

IBAM 

IBAM 

ISAM 

ISAM 

ISAM 

IBAM 

ISAM 

IBAM 

IBAM 

ISAM 

ISAM 


ISAM 

ISAM 


IBAM 

ISAM 

ISAM 

ISAM 

ISAM 

ISAM 


ISAM 

IBAM 

IBAM 


IBAM 


IBAM 

IBAM 


IBAM 

IBAM 

IBAM 

ISAM 

IBAM 

IBAM 

IBAM 

ISAM 

IBAM 

ISAM 

ISAM 

IBAM 

ISAM 

IBAM 

ISAM 

ISAM 

IBAM 

ISAM 

ISAM 

IBAM 

IBAM 

IBAM 

IBAM 

IBAM 

IBAM 

ISAM 

IBAM 

IBAM 

IBAM 

IBAM 

IBAM 

IBAM 

ISAM 


CAFTI 1)1 *»LAO 
••Mill,  1  1  MLAO 

FORMAT*//.  !•!,  ‘UNtVAAlATI  BSNBITV  ESTIMATION  ABSULTS  Ml  * 

M  •  »•  I  .N 
•Ml  I  !•«(  I  ) 
l  CONTINUE 
»•««> 

MaMORO* 1 
IFIM.ST.I)  MI| 

MMt*M* I 
Mi  1 . /•lllflMl 
irillMT.It.tl  II  »•  II 
CALL  OUtCRlN.ON) 

>  CONTINUB 


HI  1  |»  .  S«N 
•  0  SO  J •  1  .  • 

>  CONTINUE 

COMFUTB  LITTLE  •  AMD  F0«1/|L|TTLB  0) 
MIHMt 

CALL  •TIPIIIN.tt.HM ,  OL .SFCFACI 


COMMUII  WCI8MTCD  SFACIN6S  ILITTLB  0 1 U  > I  D 
00  ••  I ■  1  .  M 

>  «K1|  |  |*OOONCtM|l«l|  .  IOONI 

CALL  WSFACBI«RS.C«»S . MO  1 ,01 .WAt .U.SISO) 


IMFUTB  F0UA1SR  TRANSFORM  OF  WfICMTBD  BFAC1 
CALL  FOMIIftlWat , Wt 1 1 .M.A.M) 


H 


I  !•« 

00  III  Ri I .MMf 
RFIIM1 

CALL  AUTOABI A.RF1  .M. ALFN. FNI , MOO V AO) 

OF AM | K I  • MB AL  I  AO 0  VAR ) 

I  LOC 1*1*1  I 
00  OO  Jal.A 
ALFNAI  I  I  I • AL AMI J ) 

I  I  •  I  I  •  1 

9  CONTI  NUB 

>  CONTI  NUB 

CALL  FAAItOVAO.M-1 .N.CAT.NVW) 

IMN««.BI.I)  00  TO  110 
LOC* ILOCIMV*) 

DO  no  |||.«*» 

ALFA  I  1  I • AL  FNA( LOC ) 

LOC*LOC*t 
9  C0NTIMU8 

I  CALL  CLALT 1  I AVAA.M-  I .  1 . ANA VAN. A  I .  I  I 

COMFUTB  UN  I  TAN  1  AT B  OtMO ITT -OUANT I  LB  AT  lOO  F 

WOIVKIO. 1SOI  t IBO 
9  F  ONMAT (/.  Ill ,  *01 80  •  *. 010.81 

NVAAAi AVAR I MTW ) 

DO  IOO  I » 1 . 1AO 
U(l|*FLOAT( I  1/101 ,| 

Olal  .0 

IFIMVA.6T  O)  01  a AA88T | U I  I  ) . N VW . A V ARM . AL »N I 
1010  I  .00.0.  )  F| aH 

FO|  UaFOFMClHIII  .  I  DON  I  /  I  0  I  •  8  1 60  I 

>  C0MTIMU8 

1F| IOLTS  00.1) 

♦CALL  OLOTNTfH.OO. 1OO.CA0T, AN  U . AM  F0.1) 

BNO 

SUOAOUT I  MB  AUT0A6IA.L00  1 .M.AL0MA.0N1  ,FNN) 


COMFUTB*  THf  COBFFICIBMT*  ALFNAI.)  AMO  FAN  OF  TNB 
AUTOAf CRESS  I  VC  8*  T I  MAT  ON  ACCORD 1M6  TO  A  RECURSIVE 

INFUT  r 


OF  0 1 MB MS  I  ON  AT  L  8 AST  M 
M  :  I M *  1  )  18  TNB  MANIMUM  OAOBA  OF  8CNBM8 

TO  OB  COMFUTBO 

l OF  I  :  0A08N  OF  8CNBM8  0BIM6  COMFUTBO  FLUB  1.  L8FI.6S.B 

ALFMA  :  VB C TOR  OF  C08FF I C I BMTS  0BF)N1N6  TNB 
AFFROI I MAT  INC  FUNCTION.  NAS  TO  OB  OIMBN* 

SI0M80  AT  LBAOT  M  AMD  OBCLAABO  C0MFLC8 
FAN  :  SCALES  TNB  AUT0AB6RBSS I  TO  88TIMATOA  TO 
I  NT B6 MATS  TO  Alt).  OBCLAABO  NBAl 
ALFMA.  FNI  ANO  FAN  ARB  USBO  ABCUA8 I  TBIT , 18 .  TNB I A 
VALUBO  AT  OUTFUT  000  OAOBA  J  AAB  U8B0  AO  INFUT 
FOR  ORDER  I J* 1  ) 


COMF1!*  A | LSF1 ) . ALFMA ILSFI I.FNI (LOFf ) .8. 
COMF1  Bit  FAN 

IS  »  LSF I •  I 

TWOF 1*1. • AT AN I  1  .•) 

F JNtCMFLHlO. ,0. ) 

FMIILS1  *  CM0LBI1..O.) 

IOILS  .80.  D  FAN  ■  C  ON JC I  A ( 1 ) ) 

DO  A  I  ■  I  . LB 

A  FJMaFJMaCONJOf Af I* 1  )*FNI  |  I  )  ) 

8  •  • F  JN/OAN 

ALFNA (LSI  •  B 

IOILS  .80.  1)  80  TO  8 

A  a  LB  -  I 

00  S  I  •  I  .N 

S  ALFNA |l)  *  ALFNA (l|  ♦  CaFNI(l) 

0  CONTINUE 
00  S  I • 1 .IS 

9  FNI II)  *  C0NU6I ALFNAI LS«1  • I  I  I 

FAN  a  FAN  •  F  «#N*  CON  JB  |  F  JN  )  /  CO  NUB  I  FAN  ) 

BNO 

SUOAOUT I MB  CLFLTIlR.N, I N I T . MAMA , MM. I  OF T ) 


SUOROOTINB  TO  OOINT  ANO  ONINTBA  FLOT  TNB  N-VBCTOA  I. 
INFUT  t 

N.l 

IN I T  :  FA | NT 80  INOBB  OF  F I AST  FA (NTBO  I 
NAME  :  A  CNAAAC  TEA  LITBAAL  CONSTANT  BITIN8 
L ADSL  FOR  ■ 

MM  t  NUMOBA  00  COLUMNS  IN  FIST  (IB.TOtt 
IOOT  :  I.S  (OOINT  ON  OAA  FLOT  | 


SAM 

SAM 


BAM 

SAM 

SAM 

SAM 

SAM 

SAM 

SAM 

SAM 

SAM 


SAM 

SAM 

SAM 

SAM 

SAM 

SAM 


SAM 

SAM 

SAM 

SAM 

SAM 


SAM 

SAM 

SAM 

SAM 

SAM 

SAM 

SAM 

SAM 

SAM 


SAM 

SAM 

SAM 


SAM 

SAM 

SAM 

SAM 

SAM 

SAM 

SAM 

SAM 

SAM 

SAM 


SAM 

SAM 

SAM 

SAM 

SAM 

SAM 

SAM 

SAM 

SAM 

SAM 

SAM 

SAM 

SAM 

SAM 

SAM 

SAM 

SAM 

SAM 

SAM 


SAM 

SAM 

SAM 

SAM 

SAM 

SAM 

SAM 

SAM 

SAM 

SAM 

SAM 

SAM 

SAM 


SAM 

BAM 


SAM 

SAM 

SAM 

SAM 


SAM 

SAM 


SAM 

SAM 
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•  AM 
SAM 


GAM 


SAM 
SAM 
•A* 
SAM 
•  AM 
SAM 
SAM 
•AM 
SAM 


SAM 
•  AM 


SAM 

BAM 

•  AM 
SAM 
•AM 
•AM 
SAM 

•  AM 

•  AM 
SAM 

•  AM 


SAM 
•  AM 


SAM 
•  AM 


•  AM 
SAM 

•  AM 
SAM 


SAM 
SAM 
SAM 
•  AM 
SAM 
SAM 


•  AM 
SAM 
ShM 
SAM 
GAM 
SAM 
SAM 
SAM 
SAM 
SAM 
SAM 
SAM 
SAM 
SAM 


SAM 

SAM 

SAM 

SAM 


SAM 
SAM 
•  AM 


SAM 

•  AM 

•  AM 

•  AM 

•  AM 
SAM 
SAM 


SAM 


10  AM 
ISAM 
HAM 
ISAM 
It  AM 
10  AM 
10  AM 
10AM 
10AM 
10  AM 
10  AM 
10AM 
10  AM 
10  AM 
10  AM 
10  AM 
10  AM 
10  AM 
10  AM 
10AM 
10AM 


10  AM 
10  AM 
10  AM 
10  AM 
10AM 
10  AM 
10  AM 
10  AM 
10  AM 
10AM 
10AM 
10  AM 
10  AM 
10  AM 
ISAM 
I  0  AM 
10AM 
10  AM 
10  AM 


10  AM 
10  AM 


10  AM 
10  AM 
10  AM 
10AM 
10  AM 
10AM 
10AM 
10AM 
10  AM 
10AM 


10  AM 
10AM 
10  AM 
10  AM 
10  AM 
ISAM 


IUIMOCHAM  TO  COMPUTt  MAXIMUM  INTAO’V  CSTIMATCt  OP  010AM 

TNI  SIVARIAT1  OtPINOtNT  OfNSITV.  010AM 

010AM 

TM00K  PITTtO  MOOt l 0  WILL  00  SfTUNNtO  WITH  COirPICIINTO  010AM 

|«  COr.  VANIAOLS  NAMtS  (INOIC00)  IN  INO.  TNI  riAST  LOCATION  010AM 

or  coirncitNTS  and  inoiccs  rot  tmc  snd  moocl  in  monthhi,  i.t.,  oisam 

cor  IMtNTI  1  )•!  )  CONTAINS  COIPMCI0NT  NUMOIK  ONC  OF  TM|  ftCONO  OISAM 

■  CONCStlON  MOOt  L  CONABSPONO I NO  TO  INDI1  I N0 | Mf NT | 1 | ♦ I ) .  |TC .  010AM 

OISAM 

INPUT:  N.M  *  SAMPLI  tilt.  UN  I  VAN] ATI  MAXIMUM  OOOIN  | M« • S  010AM 

UStO  rOA  ‘0  t  VAN  I  ATI  MAX  0NOBN)  010AM 

NMI  -  COVAN  I ANCt  MATNIX  010AM 

IONO  •  VICTON  Or  ON0INIO  INOICtt  rON  StOUINTlAL  NtONtSSION  010AM 

OISAM 

AUXILLIANV:  NVAN.NVAA.0tfT  -  VtCTOAf  AND  MATNIX  OISAM 

PNOM  NOUTINI  CtONtC  010AM 

OISAM 

OUTPUT:  COP. INO  •  0tt  A0OVK  010AM 

OISAM 

SUONOUTINtS  CALLIO:  C0ONI6 . CSWttN . CLNtT t  010AM 

OISAM 

•  •••••••••••••••■•a* . OISAM 

010AM 

0  I Mt  NS  ION  MINT  I  3 )  ,  INDffT)  010AM 

0  I Mt  NS  ION  10*0(40  I .NVANI4S) ,  INOVf  1331)  .*VA0|4S)  0I0AM 

COMNLKX  PHI (00.001 .COP |#T| .O0ST I  13901  OISAM 

MM«M»M  OISAM 

MMI iMM* 1  OISAM 

L  a  MOO  I M , 9  I  OISAM 

ML • I M*  L  ) / 9  OISAM 

TWOPlat .••ATAN( I .01  01 OAM 

OISAM 

CALL  NOUTINI  CS0NI8  TO  PtNPONM  StOUINTlAL  NtCNtSOION  ON  ONI  010AM 

OISAM 

CALL  CSONtCIPMl .SO.  ANA,  IONO. 00IT,  INOV.  1990, AVAN, NVAN.NIVN)  010AM 

CALL  CLPLT « | AVAN. Nl VN. I . 4MNVAN, A1 , I )  OISAM 

OISAM 

PLACt  COt PP I C I INTO  IN  COP  PON  tACM  ONOIN  OISAM 

OISAM 

OISAM 

LOC«l  OISAM 

00  I SO  It  |.3  OISAM 

N a MINT ( I  |  010AM 

IPtX. 10.01  CO  TO  ISO  OISAM 

R I ■ HP AN |K)  OISAM 

00  I TO  NR ■ t , A  OISAM 

INDUOC)  a  INOV  INI  I  OISAM 

COPUOC  I  aOISTlRI  )  OISAM 

RiaRlPl  OISAM 

LOCaLOCPt  OISAM 

TO  CSNTINUI  OISAM 

SO  CONTINUI  OISAM 

NITUNN  OISAM 

INO  OISAM 

SUONOUTINt  CSONIOIA.NOIM.NIV,  1OND.0I0T.  1 N 0 V , MO  I M . N V A N . N V AN . N I VN )  OISAM 

. . OISAM 

OISAM 

SUOPNOONAM  TO  PtNPONM  SIOUtNTIAL  NltNtSSION  US  I  NO  CSVANlANCt  OISAM 

ON  CONNt L AT  I  ON  MATS  1 1  A | N I  V* I . N I  V ♦ I | .  OISAM 

OISAM 

INPUT:  A  -  CSVANlANCt  MATNII  (COMPLtt)  OISAM 

NOIM  •  S0«  0  I Mt NS  I  ON  OP  A  IN  CALL  INC  PNOONAM  OISAM 

NIV  -  NUM0IN  0P  INDIPtNOtNT  VANIAOLIS  OISAM 

IONO  -  INTIOtN  VICTON  CONTAININO  INOICtt  OP  VANIAOLIS  OISAM 

IN  TNI  ONOIN  TNIV  ANt  TO  It  INTINIO  INTO  TNI  MOOt L  OISAM 

MO  I M  •  0 1 MINS  I ON  OP  OIST  IN  CALL INO  PNOONAM  OISAM 

OISAM 

OUTPUT;  A  •  SMIPT  CSVANlANCt  MATNII  OISAM 

OIST.  INOV  •  Vt C TONS  OP  SUOStT  I NPONMAT I ON  OISAM 

ISAM 


CALL  MINMASIA. I A . M.N . 0M|N . XMAS ) 

PACT  a | RMAX •NMINl/IO . 

00  SO  I  a  « .M 
IP L 40*0 

00  10  Rat, SO 

L1NI 1 |S|a IOLANR 
L  1  Nt  9  I R I • I  CLANK 

10  I  INtSlKla IOLANR 
00  T0  jal.N 

IT|M»a IP  III ||A( I .Ul-SMIN I/PACT )»  .S)« 1 

INIXTt | I T IMP • 9 1/3 

IP | ITtWT .10  1  ISO  TO  1 1 

IP  t ITtMP  .St . IS  ISO  TO  10 

1P| I ITtMP . to. 1S| .00. I ITIMP .to. 1TI ISO  T0  IT 
|P|  I ITOMP .10.  14)  .00 .  I ITIMP . t0.  IS  I  ISO  T0  IS 
L INI  I  I J) a  1 1 VMS | I Nt XT ) 

SO  TO  TO 

11  L INI  I tU|a ISVMOI 1 » 

IS  IPLASaf 
IPLA01 a  I 

LINt 1 1 J) • ISVMSI4) 

LINI3I J)a ISVMOIT) 

L I NISI J I ■ I 0 VMS 10 ) 

11  |PLA0«1 

LINt 1 1 J)* ISVMSI0I 
L INtll J) • ISVMOIT | 

CO  TO  TO 

It  IPLACat 

L I Nt I  I J I ■ I S  VMS |0) 

LINI9IJ)*! f  VMS ( 3 ) 

TO  CONTINUI 

U • 1 P  L  OAT |  I  )  •  O . I)*0I VM 
WN I Tt (0.001 IU. LINI1 
IP | I PL AO. 10. 1 )WN I Tt | 4 . 009  I L I Nt9 
IPIIPLACI  .10.1 )WR I Tt 1 1 . tO*  »  L INI  3 
40  CONTINUt 

INOtXa  (N*tl/t 
00  « 1  I  a  1 .  | HOIK 

01  X HUM I  I  I  a | 0 . • P L SAT  I  I  ) -• . 1 1*0  I VN 
NINO • INOtXa  10 
00  1 9  I  a  1. NINO 
02  L 1N1NI  I  I • ISVMOI 3) 

00  03  I  a  1 t .NINO.  10 
01  L I  NON (  I  I  a  ISVMOI 0 1 

IP  I  NINO . 81  .  1 10)60  TO  tt 
00  04  I  • NINO  .  I  lO 
04  L INSNI I ) a IOLANR 
At  WNITtl t.OOlli ININ 

WN] Tilt. 004) | X MUM | I  ) .  I  a  1 ,  INOtK ) 

NtTUNN 

f 0 1  PONMAT I  IN  .2S.P0.3.9X.  * ! * .  «s. 001 A1 .  IX) | 

002  PONMAT |  IN«  .  I 1X.B0|A1  .  1  X  I  ) 

001  PONMATC1N  . OX . LIOA 1 ) 

004  PONMAT (IN  .  191. 1 1  ISO .PS . 9)  ) 

INO 

SUONOUTINt  CPTtNTIN.M.PNl .  IONO.  I  NO . C OP . Mt NT ) 


I 


lUi 


I 


•  It . 

c 

•  ISAM 

•  it  . 

c 

•  ISAM 

•  it. 

c 

1  •••  •  VICTM  CIMMIIIHI  ItllCBt  ••  MIIMiai  It  •••«• 

•  ISAM 

in  ■ 

c 

•  MAT  TNBT  WIM  CNTCtCt  «••«  VAltlt  CAUtINt 

■  ISAM 

in 

c 

» INCUL At|T tit  IN  A  SMI T  TBO 

•  ISAM 

•si . 

c 

•  ••AM 

•ss. 

e 

•  ••AM 

•1). 

c 

•  ISAM 

•at . 

c 

SUIMMMNI  CALKS;  CtWtIP 

•  ISAM 

•ii. 

c 

•  ISAM 

•  at . 

c 

•  » 

•  ISAM 

•at. 

c 

•  ISAM 

•at- 

COMP III  A|NBJM, NOIM) .tlST|MOIM) 

•  ISAM 

•a« . 

6 | Mt M$ 1  OH  INST IMtlM) .  lt*t(ll|V|  . AvAAINtV)  , NVA*(NI V) 

■  ISAM 

•a*. 

•  ATt  T  SL  /  1  .§-••/ 

•  ISAM 

•si . 

NVilllV*  1 

•  ISAM 

•as. 

•  ISAM 

•as. 

•  ISAM 

•at. 

LOC»  1 

•  ISAM 

•at . 

LCIH 

•  ISAM 

•  at . 

•  ISAM 

•  at. 

«•  1 

•  ISAM 

•at. 

tt 

■  ISAM 

•a* . 

•  ISAM 

Ml. 

T •• T • AL ( A  1  It.  10)  )*«a«AIMAt| A{  10.  1 • I  )••* 

•  ISAM 

•  tl  . 

1  P  1  TUT  .  L I  .  T |L  1  00  TB  At 

•  ISAM 

•«t. 

CALL  CBMlBPf A.NBIM.NV. IB. IB) 

•  ISAM 

•  4S. 

•  VAN|S)«ABAMA|NV.NV|  l/TAS 

•  ISAM 

•tt . 

00  SO  AN • 1 . A 

•  ISAM 

Ml. 

KID* 10*01 AK ) 

•  ISAM 

•  tt. 

•  ISAM 

Ml  . 

•  ISAM 

•  tt  . 

LCULCIM 

•  ISAM 

•  tt  . 

at 

INOVl L0C)«A10 

•  HAM 

III. 

OCf  T  ( LOC)*-A(NV,KID) 

•  ISAM 

•  •  1  . 

IOC • L  OC ♦ 1 

•  ISAM 

tta . 

at 

COHTIMUC 

•  ISAM 

••a. 

SO  TO  AO 

•  ISAM 

•tt . 

to 

■  ISAM 

••• . 

00  SO  I'K.NITN 

•  ISAM 

•  M. 

10*01 1  I  ■ 1 0 AD (  1 ♦  1  ) 

•  ISAM 

■tv . 

•t 

COHTIMUC 

•  ISAM 

•  It. 

SO  TO  lO 

•  ISAM 

•M  ■ 

tt 

MA«t 

■  ISAM 

Ml. 

IP  1  ROUMT . Lf . HI  V )  CO  TO  lO 

■  ISAM 

•  tl  . 

•  ISAM 

•  tt. 

CNO 

•  ISAM 

•  ta . 

•  ISAM 

•  tt . 

•  • 

■  ISAM 

•  tt . 

■  ISAM 

•tt. 

•  ISAM 

•  tt. 

•  ISAM 

•  tt . 

•  ISAM 

••• . 

INPUT  : 

•  ISAM 

•tt . 

A.N.K1  . K  a 

■  ISAM 

•t  i . 

•  ISAM 

•ta . 

•  ISAM 

*ts. 

OUTPUT  : 

S  ISAM 

•tt 

A 

•  ISAM 

*t». 

•  ISAM 

•  tt . 

turnout  1 Nil  CALICO  :  MONK 

B  ISAM 

•  tt . 

■  ISAM 

•ta. 

c 

•  • 

■  ISAM 

•t*. 

•  ISAM 

Ml 

■  ISAM 

•  •1  . 

DAT*  HOUT/t/ 

•  ISAM 

•  ta. 

•  ISAM 

•ta. 

•  ISAM 

•tt . 

•  ISAM 

•  ••■ 

00  SO  A ■ A  1  .AS 

•  ISAM 

•  tt 

•  ISAM 

•  tt  . 

■  ISAM 

•  tt  . 

•  ISAM 

•  •• 

TttT  •NIAL  (  A  (  A  ,  A)  ) » • 2* A  1  MAC ( A l A . A  )  )■•! 

•  ISAM 

**t. 

IP  t  T«*T  .LT.  1  . ■ -  SB )  00  TO  •• 

■  ISAM 

•  at . 

D« 1 . /A(A.A) 

■  ISAM 

••a. 

A  |  A  ,  A  1  •  »  . 

■  ISAM 

••a. 

c 

■  ISAM 

••t . 

c 

ITH  HOW  : 

•  ISAM 

•••. 

c 

■  ISAM 

••• . 

DO  It  1 •  1  .  N 

•  ISAM 

•tt. 

•  t 

At  A.  !  ) >D"A( A.  1  ) 

•  ISAM 

»*t. 

c 

•  ISAM 

»•* . 

c 

ITH  COLUMN  : 

•  ISAM 

IttO. 

c 

■  ISAM 

IMI  . 

00  ao  J» I ,N 

•  ISAM 

itoa. 

IP! J.IO.K)  CO  TO  SO 

■  ISAM 

itta. 

At J. K) • * A< J. A)'0 

■  ISAM 

ittt . 

ao 

CON  T 1 NUC 

•  ISAM 

itta . 

c 

•  ISAM 

••at. 

e 

■  ISAM 

ittt . 

c 

•  ISAM 

tttt . 

■  ISAM 

ittt . 

•  ISAM 

itit. 

do  a#  m.N 

•  ISAM 

•  tn. 

IPI 1  to. A)  CO  TO  ao 

•  ISAM 

•  •«*. 

A(j,  I  ) • A  |  J ,  |  |PA( J , A )*A| A ,  1  )/0 

■  ISAM 

it  ia . 

ao 

CONT INUC 

•  ISAM 

ttit . 

to 

CONT INUC 

■  ISAM 

toia. 

c 

•  ISAM 

•  tit. 

c 

•  ISAM 

•tit . 

to 

CONT INUC 

•  ISAM 

•tit . 

c 

■  ISAM 

•tit. 

e 

■  ISAM 

••at. 

CO  TO  1  IO 

■  ISAM 

•  tai . 

•  • 

WA ITCINOUT. 100)  A.A1.AI 

•  ISAM 

itta. 

it# 

POAMAT |  ItC,  1  a .  IONTN  0 1  AC  OP  P*OM .  1 1 .  1 7 .  1 > . INTO .  1 K , 

•  ISAM 

•taa. 

lia.tR.ITHIC  ICAO  IN  CSWCCP) 

•  ISAM 

•tat . 

1  ft 

ACTUAN 

•  ISAM 

•taa. 

CAN 

■  ISAM 

•ta*. 

•  UOAOUTIAI  DATA  1 N 1  NT AP* .  I  L  , L*t> 

•ATAIN 

•  ISAM 

•  tat. 

DATAIN 

•  ISAM 

•tat. 

DATA 1 A 

•  ISAM 

•  ta*. 

•ATAIN 

•  ISAM 

its*. 

CAABI  :  LAB  III) . L AB ( 20 )  | 2tAt ) 

DATAIN 

•  ISAM 

»•*» . 

CAABS  :  SAMPLC  SISB  N.  POAMAT  1(1) . i | B ) . 1 1 S . 4R . SAA ) 

•ATAIN 

•  ISAM 

•  •as. 

DATAIN 

•  ISAM 

•taa. 

DATAIN 

•  ISAM 

•tat. 

OATAIN 

••>•. 

OATAIN 

•  ISAM 

•  tat. 

•IMCNBIBN  II 1 I.LIII . LA9IS0) 

OATAIN 

■  ISAM 

•tat. 

c 

OATAIN 

■  ISAM 

••as. 

AB AO  1  NT  APB , f )  IBB 

OATAIN 

■  ISAM 

•taa. 

i 

PBAMATiatAA) 

DATAIN 

■  ISAM 

tttt. 

HIAtlNTAPt .*)  N.L 

•ATAIN 

•  ISAM 

•tt  i . 

a 

POAMAT | |« , 41 ,SA4 ) 

■  ISAM 

•tta. 

Al AS  1  NT  AAC , L )  11(11.  1*1, A) 

OATAIN 

•  ISAM 

•tta. 

e 

OATAIN 

■  ISAM 

ittt . 

OATAIN 

■  ISAM 

•tt*. 

CHS 

OATAIN 

■  ISAM 

•tt*. 

•  1 

■  ISAM 

»•*». 

c 

•  • 

•  •  • 

•  ISAM 

I 


I 


»••• . 
»•*• 
'099 
1*11 
IMI 
(•II 
IHI 
*••• 
«••• 

•  097 
»••• 
ttll 
*•*0 
»OS» 
I0S3 
'093 
'•94 
I  OSS 
*044 

t»ll 
I  ••• 

im 

•  •ti 
'•72 
'•73 
'•74 
IOTI 
ItTI 

'•II 
'•74 
'••• 
'•II 
'•4  2 
'••3 
'Oil 
'••I 
'•II 
'••V 
'••• 
111! 
'••© 
'•O' 
'•It 
'••1 
»•** 
'Oil 
'••• 
'097 
'••• 
'••• 
'  'OO 
"•I 

"Ot 
'  '03 
"•I 
'  'Ol 
'  lOI 
'  '07 
'  '•• 

'  109 
'  I  'O 
"  " 
'"2 
'  MI 
'  '  '  I 
'  "9 
"  »• 

'  '  '7 
"  'I 
'"9 
'  '20 
"21 
"22 
"23 
'  '31 
"39 
'  '31 
"27 
"21 
"29 
'  '30 
"31 
'  '32 
"33 
"31 
'  '39 
'  '39 
'  '37 
'  '31 
'  '39 
'  '«• 
"I' 
"12 
'  '13 
"II 
'  '19 
'  '19 
"17 
"•9 
"49 
"90 
"9' 
"92 
"93 
"91 
"99 
"99 
'  '97 
'  '99 
"99 
'  '90 
"9* 
"93 
"93 
"94 
"99 
"99 
"93 
"99 
"99 
»  130 
"  7  • 

"73 

"73 

"74 

"79 

"79 

"77 

"79 

"79 


9:  MM'  OP  99994  9747|9tlC9 
4;  0IM9M919M  91  19947  9 

N4M«  .  M4MI  OP  9474  997.MU97  99  44447  OP  PIM99II0N  •  • 
64"  INS  PI9C44M 
I "940 :  M| 49 1  49  P94  4M4179I9 

I  UN  I  7  :  NUM994  OP  OMIT  OUTPUT  It  0991490  ON. 

1417:  •  P  04  PI497  94",  '  7M94C4P794 

10U7:  I  IP  0U4MTI1I9  70  09  tISTBD,  •  9TN94WI9I 
OUTPUT:  P4I4T99  OUTPUT  19  04  I U4 I T . 


IN 


9IM949I04  1(4) . 44MC ( 3« ) , 9UM( 4 ) . 6UMS0I4 ) 

D I  Ml 49  I  04  AIP| J) .414  ) 

9  I Ml 41  I  04  I M| 40 ( 30 ) 

D4T 4  417/04 .  .  'O.  .39/ 

4UNIT  .  | 

C  COMPUTI  i.TMt  44447  9P  9U4RTI19  91299 
IP  I  IN"  .IO.  •)  SOT 0  9 
1PUI4)  .90.  4)  SOTO  39 

9  11  •  4/4 

11  •  1  L 

12  •  11 

13  •  11 

14  •  11 

19444  ■  MOO  1 4. 4)  *  ' 

SOTO  ( 30.  "  .  '2 .  '3 )  .  11444 
"  COOT  I 4UI 
14  •  11  •  * 

60  TO  30 
'3  C0N7I4UI 

1 »  •  11  ♦  » 

14  •  II  «  ' 

60  TO  20 
'3  C0NT14UI 

12  •  It  •  1 

13  ■  11  *  I 

14  •  11  ♦  » 

30  C0NT14U9 

11  ')>!' 

1 1  3  I  • 1 (  '  I  ♦  12 
1131*113)  ♦  13 
i|4)«l(3)  ♦  14 

C  POINT  04T  4  44047  •  ONI  C01UMN  POO  I4CN  OU44TIO. 

36  WO  I  T9 | NUN I T .  '•• '  )  N4MC 

WOITSI NUNI 7 . '020)  INI40 
IP  I (OUT . 90 . O)  SOTO  39 
W01T9I4UNIT. '0021 
WO  I  T9 | NUN  I T ,  '003) 

DO  3*  I  •  '.11 

30  WOITIINUNIT.  '004)  I.  9(1).  Rill"  ♦  I)  .  9  |,l  I  2  |  ♦  I  )  ,  1 1  1  I  3  |  ♦  I) 

WO  I Tl I  NUN  IT ,  '0091 

IPll'  .67.  11)  WOITIINUNIT.  'OOII  l|l(") 

I  P  1 1 3  .ST.  11)  WOITSI  NUN"  .  '007  I  91113)1 
IPI13  .07.  11)  WOITIINUNIT. '009)  ItllS)) 

I P 1 1 4  .07.  11)  WOITIINUNIT.  1009  I  14.3(1(4)) 

36  I P I  I  4 1 T  .10.  ')  OITUON 

C  COMPUTI  440  POINT  OI9COIPTI7I  9T4TISTIC9. 

«  ■  ' 

6  •  O. 

660  •  O. 

DO  90  I  *  '.4 
91  •  0.0 
90)  •  0.0 
90  »  11  I  I 
00  40  J  ■  R.KK 
SI  •  11  ♦  HIU) 

901  •  901  ♦  9lJ)*9(«fl 

4*70111) 

9  •  9  •  9 1 
990  >  990  ♦  901 
9UM| I |  •  91 

9UM90")  ■  901 

90  CONTINUK 

1  P I IOUT  .90.0)  90  TO  9ft 

WOI  Tl  INUNIT  ,  lO'O)  I  SUM  tl).M'.4) 

WOITIINUNIT.  'O"  )  I9UM90  II).  1*1.4) 

99  3940  •  9/P104TIN) 

740  >  |990  *9*I04N)/P104T|N* 1  ) 

90  •  900  T ( 740 ) 

039  •  OP  I  NO  I  3, N.  .39) 

090  *  OP  I  NO  I  3. N.  .90) 

079  •  OP  1401  * .0.  . 79) 

010  •  979  -  #39 

■940 10  ■  IRON#  -  OOO I  /  I*.  •  019) 

9010  •  00  /  13.  •  0101 

901019  •  419919910) 

SOON  *990/71047 (Ml 

T0IM4  •  I  039*3 . *090 *07 9 ) • . 29 

#43  77  •  .  3*  OP  I  NO!  3.N,  *  .  /  J )♦.  4* 090*  .  3*  OP  1  NO  I3.N.8./3.  ) 

WO  I  79 ( NUN I  7 .  '0131 

WOITIINUNIT,  10  12  1  4 . 039 . 090 . 079 . 0 1 0 . TO  I MN . 649 T 7 
WOITIINUNIT.  1 0 ' 4 )  S904 . 1940 . 740 . SO . 1940  1 0 . 19  1 0 , 90 1 919 
WOITIINUNIT . '019) 

00  40  I  ■  '  .  3 

15  •  I NT | 41 P | I )*P 104T |N|) 

N  •  N  >3*16 

0M6M'  •  N  -  IS  -  ' 

ISP3  •  IS  ♦  3 

TOM  •  R  I  19"  "SIN- 19) 

WMN  *  91 19*' )*P194T ( 16)  «  3(4*16) • P194T (19) 

00  TO  J  •  I0P2.NMSM' 

TOM  i  TOM  *  llUl 
WMN  •  WMN  •  I ( J I 
70  CONTINUt 

TOM  «  TOM/N 

WMN  ■  WMN/Pl 04 T |N) 

WOITIINUNIT.  '0171  4lPm.WWN.TNM 
SO  CONTINUI 
099  CONTINUI 
99TU4N 

'00'  POOMATI ///T30.30A4) 

'002  P0NMN7I/T4O. *00010  ST4TI9TSC9  IN  0U4NTIN9 * /T40 . 20  I  IN* | ) 

'003  P99M0T | 72# . *  9BOUINCI • /TOO. *  Wl TO I N * /T 30 . *  9U49TI1I 

♦  * P I  097  0U49 TIN  9IC0N9  904KT90  TNIND  9UART99  POUOTM  OU40TIR* 

*  /T31.0I 'N* ) .3(33. '3| 'N* 1.39. 141  IN*) )  ) 

'004  PONMA7I7JO. 10  .  4(11.711.4)) 

'001  PONMOTI '3) 

'009  POOMATI IN* . T29.P '9. 4 ) 

'007  P99MAT | ' N* ,749, P'9. 4) 
f 009  P90MATI IN*. T9' .P19.4) 

'009  POOMATI IN* .730, 14. TT7. P'9. 4) 

'O'O  POOMATI 'NO. 730. *9UM* .93.41 '3. P19. 4)) 

'O"  P9AMATI 'NO.TSO. 'SUM  9P*/TS9. *  S0UARI9 * , 1 9 , 4 1 1 3 , P 1 i  .  4  )  | 

'0*2  P0NMAT I///T49. *  999CN I PT I  01  9TAT 1ST ICt * /T49 , 33 1 IN* ) J 

'013  POOMATI//*  SUMS TAT  * . 730,  *  9 AMP 19  * , T2I ,  *  1 0W9R * . T98 . *  UPPIR 

*  T  94 , *  I NT  OU ART l * . T9S . *  S4STWIRTNS * / *  SUMS TAT *. 730 .  *  till* 

♦  .739.*  •UARttiC* .749. *  MID  I AN  * . T9 1 . *  9UARTI19*. 

•  T94 . *  RANSI ’  , T7  4  .  *  TRIMtAN* .799. *  ttTIMATt* 


9AM 

SAM 

SAM 

SAM 
9  AM 
SAM 
SAM 


SAM 

SAM 

SAM 


SAM 


•  AM 
SAM 


SAM 

SAM 

SAM 


SAM 

SAM 

SAM 

SAM 

SAM 


•  AM 


SAM 

SAM 

SAM 


SAM 
•  AM 
SAM 
SAM 


•  AM 

SAM 
SAM 
SAM 
9  AM 


SAM 

SAM 

SAM 


SAM 

SAM 


9  AM 
SAM 
SAM 
SAM 
SAM 
SAM 
SAM 
SAM 
SAM 
SAM 
SAM 
SAM 
SAM 
SAM 
SAM 
SAM 
9  AM 
SAM 
SAM 
SAM 
SAM 
SAM 
SAM 
SAM 
SAM 
SAM 


SAM 

SAM 

SAM 

SAM 


9  AM 
•AM 


ISAM 

ISAM 

SAM 

9AM 

9  AM 


•AM 
9  AM 


SAM 
SAM 
9  AM 


ISAM 
SAM 
ISAM 
9  AM 

SAM 

•AM 
•  AM 
SAM 
SAM 
SAM 


•  .//■  CUMSTAT*  .TOO.  It .  T39  .§<Sft  .0,  Ill  I 

t*IA  FAAMAT|///‘  SUMS! AT* .TfO, ‘tVMOO/N* , *  MBAN'.TA#.*  TAB I AN  CO * , T ftt . 

•  •  ITI  IIV  .HI.'NIM  ll'.m/ltl  IIV  ll'.TII. 

«  *  ICC  0TB  10*//*  •UMIfAIMU.YIIII.I.IIlH 

•  O 1 0  P  OAMAT  |  / / / T  40 ,  *  T IIIMIC  AT  I  ON  *  .TBS,  * Ml  NO BA I  SSB  ',TM,  *  TAIMMB9’/ 

•  U*,*  POINT*. TBS. *  MOAN*. TOO.*  NCAN'/I 

t« IT  P OAMAT I  TOO, PT, 3. Tit, 9(01 I  .4.1*11 

1090  FOAMAT IT90.90A4) 

■  NO 

. . .  * 

t  SWiAOUTIHI  TO  COOVBAT  AIAI  VAAIAitO  > 

C  W« 1  CM  MAO  «  CHAOACTOA  V* FOAMAT  I  FOAM 

C  TO  •  C  MAR  AC  TO  A  At  PNAMB A  I  C  A  ABA  V  NAMB  MM1CM  10 
C  IN  A  POAMAT . 

C  INPUT  ;  HOC ACM  t  OCA AT  CM  TAP*  MWMBSA 
C  I.  1P00M 

C  OUTPUT  •  NAMCITI .NAMBISI  :  A  CMAAACT t At  CACM 

. . . . . . . . . 

0  I  MO  MO  1  ON  MAMOtal.  IPOANOI 
HOC ACM  •  tl 

«A I  TO (NSC ACM. 1P00MII 
AOWINO  MOC ACM 

10  POAMATIIAAl’ 

AOTUAM 

•NO 

OUOAOUTIMO  POAIOAIP.U.H.A.MAI 


OUOAOUTINt  TO  COMPHTB  TMt  POUAIfA  TAAMOPOAM 

PM  I  I • I  OP  A  MONO  ITT  OOPIMIO  AM  10,1)  POA  V«0.1.....M 

INPUT  : 

P.U.M  !  VOCTOAt  OP  1 1 MOT M  M  COHTAIMtMC  PtUt.U 
•M  MAI I  MUM  VAtUI  OP  *  POA  MMtCN  PMIlV)  li  COMPUTED 
OUTPUT  ;  A  ;  COMP 1 1  * • VAt WOO  VtCTOA  COMTAIMIMC  TNI  PNl'O 
t UO POUT  1 000  C At  tCO  :  MONO 


OIMONSIOM  PIOI.UIO) 

COMPtA*  A | MA |.l 
T WAP  1*0 .  • AT AMI  I .  I 
PM«  P  t OAT | M ) 

00  DO  I M*  t  ,  MA 
P IM* IM- t 

A  I  IM|  aCMPtMO.  } 

00  to  IH.a 

t  » CMP 1 1 ( O ,  .TWOPl*P|««Ui I  I | 
lO  At  I M ) ■ A I  I M | *P I  I ) a  C • *P  t t ) 

A l IM) • A |  IM) /P  tOAT |MI 

to  comtimuo 

A  I  t ) aCMPtAl I .  .0.  I 
AOTUOM 
•  NO 

PUMCT10M  POPNCIB.IOOM) 


C  AOUTINC  TO  COMPUTC  TNI  VAAIOU9  DO  NO ITT  *  OUAMT I Lf  FUNCTIONS 
C  INPUT: 

C  «  •  VAL  Ml  AT  AN I  CM  TMt  FUNCTION  10  TO  01  COMPUTIO 

C  IOOM  *  INOICATOA  POA  INI  OlltMT  FUNCTION. 

C  MUST  OB  IN  TMC  KICtUOITt  A AMOS  1*1t 

COMMON  /PAAM/  OtTAP . ACT AW 
60  T0M.3.9.4.S.B.T.0.9.1O.11). IOON 
C  COMPUTC  TUB  NOAMAt 
1  CONTINUE 

IP(  *  .tT.  OOU  00  TO  lOI 

|p|  a  cr.  .tool  so  to  «oi 

PI •• . • A  TAN |  I  . O) 

CAtt  MON A IB | B, T . ISA I 

POOOaCaPI •O.B*T«»B|/IOOAT|a.«PI  |  | 

POP  NC  •  POOO 

»OI  poooao. 

PC PNC  •  POOO 

C  COMPUTA  T  MB  B  BPONBNT I At 
I  CONTI  NUB 
POOO«  »  . ♦* 

POPNC  •  POOO 
BO  TO  *0 

C  COMPUTB  TNC  tOSIOTIC 


CONTIMUB 
POPNC  •  *•!». 


A  CONTI NUB 

|P|B  IB-  .B)  POPNC  •  K 
IP  I  a  .  CT.  . B |  POPNC  •  I. 

C  COMPUTB  TNB  UNIPOOM  ABCIPOOC 
O  CONTINUB 

POPNC  i  II.  •  •)••> 

C  COMPUTC  TNB  CAUCNT 
B  CONTINUB 

PI  ■  4 . • ATAN | I . O) 

POPNC  •  OtN|P|*B)"S/PI 

C  COMPUTB  TNB  BBTOBMB  VAtUB 
T  CONTINUB 

IP  I  a  .BO.  I . I  BO  TO  109 
POPNC  •  I A  •  I.MAtOBfl. 

POPNC  •  O. 

C  COMPUTB  TNB  tOB  NOAMAt 


•  CONTINUB 

IPII  .BO.  1.1  II  TO  f 09 
IP  IB  BO.  O.)  SB  TO  103 
CAtt  MONO  IS | B , T , ISA ) 

PI  •  4 . • AT AN ( 1.0) 

POPNC  •  | •O.B4T*43“T| /BOAT  I  9 . *P I  I 
POPOC  ■  IIP | POPNC I 

103  POPNC  •  0.0 

C  COMPUTB  TNB  PAABTB 

•  CONTINUB 

POPNC  •  |I.*B|«4|1.«BBTAP)/0BTAP 

BB  TO  00 

f  COMPUTB  TNB  •CIOULt 
•O  CBNTINMB 

IPII  BO.  1.1  BO  TO  100 

POPNC  •  1 1  .  •S|»(-AtBSf  1 .  *B|  |oo|  1  . -BBT AM) /BIT AM 

BB  TO  OB 
to#  CONTINUB 
POPNC  •  O. 

BO  TO  OB 

C  COMPUTB  TNB  NAtP  tOSIOTIC 
tl  CBNTtNUB 

POPNC • I . •■••I 


I»1l. 

99  417 1)4  4 

OIOAM 

1**1. 

IN* 

014AM 

*1*4 

9  ISAM 

till. 

c« 

914AM 

n»* . 

c 

0  14  AM 

»**t. 

e 

014AM 

tsi*. 

c 

R  ISAM 

*»**. 

c 

0  1  SAM 

»s»o 

c 

INPUT:  u  •  VICTOR  or  V4LUSS  AT  WHICH  7  (VALUATOR 

ISSI  . 

c 

914AM 

tat*. 

c 

OIOAM 

isaa 

e 

0  19AM 

tsi« 

c 

014AM 

*i*i 

c 

0  1  SAM 

mi 

c 

9ISAM 

m* . 

c 

014  AM 

**»*. 

c 

OIOAM 

mi . 

c 

014AM 

m«. 

c 

1  ISAM 

**** 

0  1  9  AM 

1SSS  . 

9IM4N910N  U|N),V|N),I(M),P|M) 

019AM 

tsas. 

IP(N.IO.M)  49  TO  *04 

OIOAM 

«**«. 

ll«* 

014AM 

tut. 

90  40  I • 1 .M 

014AM 

mi 

(O  IP  III  I |-tt| II | 1*0.40.40 

II  SAM 

*nv 

lO  IMIl.ll.il  II  TO  11 

014AM 

tut 

rm«vi(}«iviii-vm)*iim-uiii)/iui2)*uiii) 

014AM 

11*9  . 

49  TO  40 

014AM 

**«•. 

30  ril|*V|ll~l|«IV<ll)»V|ll-1)t'mi)-Utll-*H/|M(lll*U(ll-1|) 

014AM 

1*4  *  . 

CO  TO  40 

014AM 

*14*  . 

•  o  rt  i  I-VI 1 1 1 

01  SAM 

1141 

CO  TO  40 

014AM 

1144 

90  !  1  »l *•* 

01  SAM 

**«*. 

ir l 1 1  . IT  .N)  40  TO  to 

014AM 

1*44  . 

1  1  »N 

OIOAM 

*147 

CO  TO  10 

OIOAM 

1144 

40  CONTINUC 

01  SAM 

1144  . 

OO  4CTU0N 

0  ISAM* 

1140  . 

•  NO 

019AM 

*14* 

9U0R0UT INK  1 C 009 A | 4 . 1 F ORM . N4Mf | 

014AM 

*143  . 

c 

OISAM 

1141  . 

c 

S  UO A  OUT  INI  TO  C0NV90T  INT94KN  VARIAOLI  K 

0  ISAM 

1144  . 

c 

WHICH  H49  4  CNANAC  TIN  1  •  4  All  MAT  1  P  ONM 

OIOAM 

1144  . 

c 

TO  4  CHANACT4N  ALPHAMERIC  ANNAT  NAMt  WHICH  19 

OISAM 

1144. 

c 

IN  A*  FORMAT  . 

0IS4M 

*147  . 

c 

INPUT  :  N4CNCH  i  4CRATCH  TAPI  NUM9CN 

014AM 

*144  . 

c 

K.  1 P  ONM 

OISAM 

*144 

c 

OIOAM 

1140 

OISAM 

*141 

DIMINIION  NAMKI2 ) . IPONMI 1) 

OISAM 

*142  . 

NSCOCHi 1 1 

OIOAM 

*141  . 

NIWINO  N4CNCN 

OIOAM 

*144  . 

WAIT* IN4CNCH. 1  FORM) A 

OIOAM 

(144  . 

RCWINO  N4CACN 

OIOAM 

(144  . 

ACAOI NIC ACM , (OINAMK 

OIOAM 

(147  . 

10  FORMAT I3A4) 

OIOAM 

1144  . 

■1 TUAN 

OIOAM 

*149  . 

■  NO 

OIOAM 

*174 

9UOAOUTINC  KINO AL IN. T AUA . T AUO . 4 0M« A . NC . NO . N 1  NO . N09P . NPA1 A4 ) 

OIOAM 

*17*  . 

c 

OISAM 

*172 

OISAM 

1371. 

c 

OISAM 

1374 

c 

OISAM 

*374 

c 

RANKS) .  AND  SOMRR'S  0. 

OtSAM 

*374. 

c 

OISAM 

1177 

c 

INPUT:  RANKS  •  TNI  VICTOR  OF  RANKS  OP  TNI  IN0IPIN9IKT  VANIAOLI 

OISAM 

*174  . 

c 

OIOAM 

*174 

c 

OIOAM 

1140 

c 

OISAM 

*141  . 

c 

OUTPUT:  TAUA  -  KINOAIL‘9  TAU- A  A94UMIN4  NO  TIK*. 

OIOAM 

*143. 

c 

T  AUO  •  KIN0AL1  '  4  T  All  •  *  POO  Till  RANK*. 

OISAM 

(141  . 

c 

••Ml A  -  SOMRR’S  0  PON  TICO  RANKS  IN  TNI  OIPINOINT 

OISAM 

(144  . 

c 

VAR  1 A4ll  . 

OISAM 

(144  . 

c 

NC  *  NUMOIR  OF  CONCORDANT  PAIR* 

OIOAM 

(144. 

c 

NO  ~  NUMOIR  OP  0 1 9C0R0ANT  PAIRS 

OIOAM 

*147 

c 

NINO  •  NUMOIR  OP  PAIRS  WITH  Till  IN  TN*  IN9IPIN0INT 

OISAM 

1144 

c 

VAAIAOif  OUT  NOT  IN  TNI  OIMNllNT  VARIAOLI 

OISAM 

1149 

c 

•01 P  -  NUMOIR  OP  PAIRS  WITH  TII*  IN  TNI  OIPINOINT 

OISAM 

1140 

c 

VARIAOLI  OUT  NOT  IN  T  HI  INOlPINOINT  VARIAOLI 

OtSAM 

*14* 

c 

NPAIR4  •  NUMOIR  OP  PAIRS  OP  0 1  VAR  1  ATI  OOSIRVATIONS 

OIOAM 

(142. 

c 

OISAM 

1141. 

c 

OIOAM 

1144. 

c 

OISAM 

*194. 

OISAM 

*194 

c 

OISAM 

1197 

COMMON  /OATAR/  1 1 900 |. T 1 400 ). RANK 1 | 400 ). RANKV | 900 1 

OIOAM 

*399  . 

0 1  MINS  1  ON  T1I4O0.3I 

OISAM 

*199  . 

00  lO  1 • 1 .N 

OIOAM 

(400 

T  1  I  1  .  1 | *040X11  |  ) 

OISAM 

(40*  . 

Till  .*| iRANKTl  |  | 

OIOAM 

*401  . 

*•  CONTINUI 

OIOAM 

*401. 

c 

OIOAM 

1404  . 

c 

OISAM 

*409  . 

c 

OIOAM 

1404. 

NC  « 0 

OIOAM 

1407. 

NO  >0 

OIOAM 

*404. 

N I N0*O 

OISAM 

1404  . 

NOIP  *0 

OIOAM 

14  10. 

1  1  1 

OIOAM 

14  1*. 

lO  K  ■  1  ♦  1 

OIOAM 

*41*. 

OISAM 

*411  . 

c 

OtSAM 

1414  . 

c 

CHICK  IP  RANKS  AIR  IOWAL 

OIOAM 

14*9  . 

c 

OISAM 

(414. 

c 

T4LCRANCI  POR  RIAL  AR 1 T NM| TIC:  AS SUM*  MO  IP  AOS | R ) . LI . * . 1  - 4 

OIOAM 

*4*7. 

c 

OIOAM 

14  19  . 

TIS  T 1 • AOS  IT*  II  .  1  I -T 1  IK.  1  I  ) 

OISAM 

*419  . 

1 P 1  TROT  1 .07 ,4.1*41  44  TO  40 

OIOAM 

1410. 

T IS  T  2 • 409 1 7 1 | |.*|-71|R,S)t 

OIOAM 

*411  . 

IPI TROT! . ST .4.1-91  N!N9*NIK0*I 

OIOAM 

*41* 

MM  1 

OIOAM 

*41*. 

90  TO  10 

OtSAM 

*4*4. 

c 

OIOAM 

14*9. 

c 

COUNT  NUMOIR  OP  CONCORDANT  AN9  DISCORDANT  PAIR* 

OISAM 

1414. 

c 

OISAM 

*4*7. 

40  00  40  UAft.N 

OISAM 

•  4*9. 

IP|TI|U,*|-TI|  I  ,  *  )  ,4T  .9  ,R-4|  NONC«1 

OISAM 

14*9. 

IPIT*|J,*)*TI| I.RI.LT. *4,1-4)  N0*N0«1 

OIOAM 

14*0. 

40  COOT IRtfl 

OISAM 

(4*1  . 

40  1  •  1  ♦  * 

OtSAM 

14**  . 

OIOAM 

•  •**  . 

c 

OIOAM 

14*4. 

c 

SWITCH  COLUMNS  OP  T 1 

OISAM 

14*9. 

* 

OISAM 

14*4. 

00  70  1 •  1  .  N 

OIOAM 

*4*7 . 

T 1 f I , » 1 •RANK I f|| 

OIOAM 

14*4. 

Tllt.1 I • RANKV III 

OIOAM 

14*9. 

70  CONTINUI 

OISAM 

1440. 

e 

OISAM 

1441  . 

c 

OISAM 

•  44*. 

c 

OISAM 

144*. 

OIOAM 

im 

*  ■  !• 1/ . 33 • t . 

•  l#AM 

t*tt. 

*  •  33 

■  i  a  am 

uta 

IM  •  1 

•  i  a  am 

Hi*. 

i  •  1 

•  ISAM 

in*. 

•  HAM 

t«»i . 

*1  •  1  . 

•  ISAM 

tut. 

I*  •  • 

019AM 

mi. 

019AM 

HM 

it  «•  •  it  •  iiiM.iiil  •  at 

OIOAM 

«•••. 

•  a*  •  »  •  1 . -*3 

#  t  a  am 

till  . 

*3  • 

OIOAM 

mi. 

a*  *  a# 

•  ISAM 

tin . 

1*3  •  *  #3  •  1 

OIOAM 

mi . 

null  .it.  h!  ••  t#  i# 

■  ISAM 

tin 

•  •  T#  1 3# . 3#  1  . i 

*  19AM 

tin 

3#  V  ■  «  ■  I#  •  *I#M* 

019AM 

mi 

••  »•  4# 

a  ISAM 

tin 

sa  •  •  #••? i •**•#( i . •*•!) i 

019AM 

HO* 

nil  .u  •*!•  )  •#  T0  3* 

0  19AM 

till 

m  •  M 1 • **M3 

0  1  BAM 

till 

IM  •  3* 

#  i  a  am 

tin. 

L  •  S 

•  IIAR 

tin 

N  •  II 

•  ISAM 

mi . 

•  IBAM 

!•••. 

3#  w  •  N3  •  a  ♦  n* 

•  HAM 

■ 

IM  •  41 

•  HAM 

»••*. 

II  ■  3* 

•  IBAM 

(••1  . 

1  ■  3 

•  HAM 

till  . 

•  HAM 

HU 

019AM 

till  . 

«•  V  •  *!•* 

•  HAM 

tin . 

■  HAM 

mi . 

a  i  a  am 

i 

IB*  •  13* 

•  HAM 

tin 

WIIITBI*.  HI  1 

019AM 

i# 1 1 

T 1  OH  IH 

///I 

•  HAM 

i*ia . 

•  HAM 

1*13  . 

•  HO 

■  HAM 

1*11  . 

MIN 

■  HAM 

1*1*  . 

•  ••  • 

MIN 

019AM 

1*1*. 

C 

MIN 

■  HAM 

1*17  . 

c 

TNI 

MIN 

•  IBAM 

1*1* 

c 

or 

MIH 

0  19AM 

1*1*  . 

c 

MIN 

019AM 

1*30  . 

c 

MIN 

■  HAM 

1*31 

c 

1MUT  i 

M|M 

•  HAM 

1*33  . 

c 

N.  XI  1  1 . ■  1  *  1 

MIN 

•  HAM 

1*33  . 

c 

Ml  H 

•  HAM 

1*34  . 

c 

OUTPUT  : 

MIN 

0  HAM 

1*3* 

c 

IMIN .  1  HO 

MIN 

•  19AM 

1*3* 

c 

MIH 

•  HAM 

1*37 

c 

M|N 

■  HAM 

1*3*  . 

c 

MIN 

•  HAM 

1*3*. 

c 

■  •»  ■ 

MIN 

•  HAM 

1*30. 

c 

MIN 

■  HAM 

1*31  . 

OIMBH9ION  ■  (  N  ) 

Ml  N 

019AM 

1*33  . 

c 

Ml  H 

•  HAM 

1*33. 

XMIN-Xt 1 1 

MIH 

■  HAM 

1*3*  . 

IHO-1 

MIN 

•  HAM 

1*3* 

IMH.IO.M  HI  TURK 

MIH 

■  HAM 

1*3*. 

c 

MIN 

•  HAM 

1*17. 

DO  1  1 ■ 3  .  H 

MIN 

■  HAM 

1*3*. 

IM1NII.01.1III]  1  HO  »  1 

MIH 

•  HAM 

1*11 

1  KM 1 H 1 1 1  1  HO  1 

MIN 

■  HAM 

!••• 

c 

MIN 

•  HAM 

I*«1  . 

■•TORN 

MIN 

•  IBAM 

1*«3 

■  HO 

MIH 

•  HAM 

1  **3  . 

4UBH0UT  INI  M  |  HMI 1  A  ,  U.M.  H.  IMIN  .  IMM  | 

•  HAM 

1*«« 

c 

•  •• 

•  ISAM 

t*«*  . 

c 

•  HAM 

1*«*  . 

c 

ROUTIHI  70  MHO  TNI  MIHIMUM  AHO  MAXIMUM 

Or  A  MATRIX 

■  IBAM 

•  *«7  . 

c 

■  HAM 

1*«* 

c 

INPUT;  H.M.AIH.M) . !*•«•«  OIMIHSIOH  OP  A 

■  ISAM 

•  *••  . 

c 

■  HAM 

«••• 

c 

OUTPUT:  XMIH.XMAR 

■  ISAM 

«**1 

c 

■  HAM 

1**3  . 

c 

•  IBAM 

t*«3. 

c 

•  HAM 

1*M 

c 

•  • 

•  HAM 

iff* 

c 

■  HAM 

1*»*  . 

0 1  MINI  1  OH  *|  1*. HI 

•  ISAM 

tin. 

*M 1 H ■ A | 1.1) 

•  ISAM 

1*3* 

■MA* ■ A | 1.11 

•  HAM 

i**t 

00  lO  1*1. M 

■  (•AM 

I*** . 

00  10  J>t.H 

•  HAM 

i**i . 

IPIAII.JI.ST.  IMA  ■  )  RMA  X  •  A  |  I  ,  tl  } 

•  HAM 

t**3. 

IPIAII.JI.IT.XMIHIXMIHNAU.J) 

SHAM 

1**3. 

lO  COHTIHUI 

•  ISAM 

t**« 

HI TUHH 

•  ISAM 

1***  . 

■  HO 

•  ISAM 

1***  . 

•  ISAM 

tin . 

c 

•  HAM 

i*ii 

c 

SHAM 

»•** . 

c 

•  HAM 

tui. 

c 

SUBAODT 1 #■  TO  S##t  TNI  HX3 -MATH  1 X  I  PROM  BMAl IBST  T# 

•  HAM 

•  it*! 

e 

•  HAM 

1*73. 

c 

•  HAM 

1*7*. 

c 

JViV  1»**1.  AHO  I#  POUHO  IH  KIRCH.  *.(1*71) 

•  HAM 

1*7* 

t 

7#RRi 

•  19AM 

1*7*  . 

c 

MOL  T  .  ■  INBMANT  .  AHO  tHIRttOH.  P.  3*7. 

■  HAM 

1*77  . 

c 

■  HAM 

1*7* 

c  • 

■  HAM 

1*7*  . 

e 

■  HAM 

11*0 

■  HAM 

1**  1  . 

na 

•  HAM 

1*9  3  . 

1  L  *1/3 

•  ISAM 

1  **3  . 

IP|L 13.3.3 

•  ISAM 

f  •••  . 

3  HI TUHH 

■  HAM 

«••*. 

3  M* H- i 

•  ISAM 

•  •••  . 

00  0  1 #  1  ,M 

•  ISAM 

tin. 

J*  1 

■  HAM 

1999  . 

4  jjtj«i 

•  HAM 

«•••  . 

trial j. il-xi jj. i III.#.# 

•  HAM 

till. 

c 

•  HAM 

lilt. 

c 

•  HAM 

1903. 

c 

•  HAM 

till. 

c 

trial  j.  i  i-ai  jj.  «i  n.a.a 

•  HAM 

till. 

c 

•  HAM 

mi. 

c 

■  HAM 

If**. 

« 

■  HAM 

nn. 

•  a i *  a ( j . «) 

•  1  •  fes? 

mi . 

8S*B| J. f  ) 

■  HAM 

i  it* . 

au,  i  MU  jj.  i) 

■  HAM 

Itll. 

au.SMBl  JJ.tl 

■  HAM 

It«l  . 

■  I  jj.  i  Mil 

■  HAM 

IMS. 

ai  jj.  tMSS 

•  ISAM 

IMS. 

J*J*t 

•  HAM 

It#*. 

in  j?*.*.« 

■  HAM 

It##. 

•  eORTIONI 

•  HAM 

It##. 

ft  Tl  « 

■  HAM 

It#t, 

BN# 

•  HAM 

•  tM  . 

•  10AM 

. 

•  •• 

•  1  RAM 

•  *•#. 

•  ISAM 

ITH  . 

0  1 OAM 

»»l». 

C  IN#V?  ; 

•  10AM 

•  TlJ. 

C  M.mail),.  ...•VIIIMI  1  OTAMOABOI  *>•  *■* 

VAR 

01  OAM 

111!  . 

C  VRR  •»••■»  1  THRU  M, 

•  1  OAM 

till 

C  «  t  OAMOtfl  IIH 

0  1  OAM 

*T»*. 

•  HIM 

•  TIT. 

c  «•••  t  •itasMinci  •*•■■ 

R  1  OAM 

»T»0. 

C  CMIll . CAT|M| 

■  I  OAM 

It**. 

C  OVOAOUTIMM  UUII  i  Ml# 

OIOAM 

IT*#. 

•  •• 

OIOAM 

IT** 

•iMCMitan  rvarimi .catim) 

•  I  OAM 

«TI*. 

ra»*w*atir) 

•  1  OAM 

•  T»». 

■  •  1  HI.M 

•  ISAM 

ITT*. 

Ci*. 

•  10AM 

ITS*. 

••  *  JO,! 

OIOAM 

IT**. 

*  C«CMt  (.•luutu  1/*#  1  t/RVAA(  JM 

•  ISAM 

1  T  JT  . 

C  •€/*# 

■  ISAM 

IT**. 

1  C*T  (M*C*I|I  .•IlLliTIII/HlI/mil  |  |  ) 

OIOAM 

IT**  . 

C*li  MIHlCAT.M. e#TMI . N*RO  ) 

OIOAM 

ITS#. 

1 V 1  CAT MI  .IT .*1 . *  * . /AN  1  HOAR •* 

•  10AM 

*  T*  1  . 

••TUAN 

*  10  AM 

ITS*. 

■  #0 

01  SAM 

ITS*. 

OIOAM 

ITS*. 

c 

OtOAM 

ITS*. 

■  10AM 

ITS*. 

c 

0  1  RAM 

ITST. 

•  IRAN 

ITS*  . 

C  COI’MCICMT. 

0  1  RAM 

ITS*. 

c 

OIRAM 

IT##. 

C  IMAUT.  I  -  INC  VICTOR  IF  VRIUOt  RF  TMf  1 MO* Ff  NOINT  V*N!**it 

it  RAM 

IT**. 

OIOAM 

IT** 

01  SAM 

IT**. 

OISAM 

IT** 

C  OUTFUT:  ■  •  FRAABRN’S  AAO*UCT  MOMINT  C*ARf l  AT  1  ON 

COCAAICICNT 

OIOAM 

»T*§. 

C 

OISAM 

IT**. 

•  ISAM 

*  T*T  . 

c 

•  ISAM 

IT*#. 

•  ISAM 

IT**  . 

c 

•  ISAM 

IT*#. 

•  ISAM 

*  T*  I  . 

c 

■  ISAM 

IT**. 

SUM! • • . 

•  ISAM 

IT**. 

SUNK#. 

•  10AM 

IT**. 

■VMM V •# . 

•  10AM 

IT**. 

•  IMIll*. 

•  ISAM 

IT**. 

•UMVViO. 

•  ISAM 

lt»T. 

c 

•  ISAM 

IT**. 

C  COMAUTI  RUMS  #A  SOVARtS  AND  CNOCf  AAOOUCTS 

OIOAM 

IT**  . 

C 

•  ISAM 

IT*#. 

00  I#  1  *  1 , R 

•  ISAM 

1  T*  1  . 

BUMI*OUMI«l|  I ) 

•  ISAM 

IT**. 

•  ISAM 

1  T  *S  • 

RUMS T • OHM! 1*1(1  1*7(11 

•  ISAM 

IT**. 

RVM11 |RUM1M«R|  1  ) • Ml  I  ) 

•  ISAM 

IT** 

R  4JMV  V  •  RUNT  V  ♦  V  I  I  )• V  (  I  ) 

•  ISAM 

IT**  . 

I#  CONTI  Hill 

•  ISAM 

IT*7. 

C 

OISAM 

IT**. 

•  ISAM 

IT**. 

c 

•  ISAM 

ITT#. 

B  *  V  • RUM1 T  - RUMN*B WMV/Fl OAT ( R ) 

•  ISAM 

ITT!  . 

R ■ ■ *  RUHR ■ *  R  UNI ■ I UMM /FtOAT(N) 

•  ISAM 

ITT* 

R  V  T • RWMVT  *  RUMV  *  RUMV / A  L  OAT | N ) 

•  ISAM 

ITT*. 

•  ISAM 

ITT*. 

RITURM 

•  ISAM 

ITT*. 

•  NO 

•  ISAM 

ITT*. 

. VMIR. VMAR. VIRC) 

•  ISAM 

IT7T  . 

RC) 

•  ISAM 

ITT*. 

MCMI  •  RC  •  1 

•  ISAM 

ITT*. 

ARC  •  1 .  /  At RAT | MCMI ) 

•  ISAM 

IT*#. 

RRA1  •  RR  •  • 

•  ISAM 

17*1 

c 

•  ISAM 

IT** 

00  10  ft  ■  t.RL 

•  ISAM 

IT**. 

IA  |  RTRTdtl  *  f  MO  lit)  1  30.30.*# 

•  ISAM 

IT»«  . 

3#  RTRTtltl  •  RMRStlSTRTl It  1 .VM1R) 

•  ISAM 

•  T#». 

*R#(ltl  •  AMIRIICRDIILl  .VMAR) 

OISAM 

IT**. 

OOTO  *# 

■  ISAM 

IT*T. 

••  *TOT (It)  •  AMIRI (OTRTI  It  1 . VMAR) 

•  ISAM 

•  T**  . 

CAR  (It)  l  A  MAR  |  UNO  tit)  .  TM1R) 

■  ISAM 

IT**. 

0#  OR  •  t  *R0 III)  •  STRT(ltl)  •  ARC  •  VIRC 

•  ISAM 

IT*#. 

*  >  l*TRT| It)  -  VMIN  )  •  VIRC  ♦  1.* 

•  ISAM 

1  T*  1  . 

INO  ■  INTI  01 

•  ISAM 

IT**  . 

•  ISAM 

IT**. 

0  •  •  ♦  OR 

•  ISAM 

IT**. 

IA  |  1  NT  I R I  .•#.  INO  I  OOTO  •• 

•  ISAM 

IT**. 

INO  •  INTIO) 

•  ISAM 

IT**. 

AANRAVlNRAI - INO. 1  I  •  CHAAt (It) 

•  ISAM 

1  T*T . 

OO  CANT  I AM* 

•  ISAM 

IT**. 

OOTO  I# 

•  ISAM 

IT**. 

SO  IA  |  Of  NT  (Hi  .tl.  VMIN  OR.  *  TNT  III)  .  Bt  .  VMAR)  OOTO  I# 

•  ISAM 

I*##. 

INO  »  (OTRTI It)  •  VMIN)  •  VIRC  ♦  1.0 

•  ISAM 

t##l  . 

00  TO  1  >  3.NCMI 

•  ISAM 

!••*. 

FANNAT INAF1 • IN*. ■ |  •  CHARI (It) 

•  ISAM 

TO  CRRT 1 AMI 

•  ISAM 

1  ••*  . 

1#  CRRT 1 AM* 

•  ISAM 

I*#*. 

•  ISAM 

•  •••  . 

«A0 

1*#T. 

OMRROUTIRC  At#TR V ll.V.N. CART ,RAMR. ARMY,  1  OAT  ) 

•  ISAM 

I##*. 

•  ISRM 

I##*. 

C 

■  ISAM 

«•!#. 

V  AS  A 

•  ISAM 

»#»  1  . 

•  ISAM 

tats. 

c 

OISAM 

ISIS. 

c  IMAttT  :  A.I.V  -  X  10  ORRIRIO  OR  1  ROUT  ARR  V(I>>V|I(I)) 

•  ISAM 

i#i«. 

*  C##T  •  lITIRAt  CRROTRRT  AOR  T 1 T t •  #A  At#1 

•  ISAM 

nit. 

C  RAMI.RAMT  *  •  CHARACTER  tITIRAt  COROT ART| 

•  ISAM 

tin. 

C  tAOItO  AOR  R  ARR  V 

•  ISAM 

•  •IT. 

C  IRAT  1  T.s  IARIRT  RR  RRR  At #T ) 

•  ISAM 

•  ••*. 

t 

•  ISAM 

lift. 

C  OttORRMT l#«0  CAtlRR  t  ATIRA . MAI . Ml R 

•  ISAM 

•  •9#. 

C 

•  ISAM 

!•«•  . 

•  ISAM 

•  •99. 

c 

■  ISAM 

•  •99. 

1) 

■  ISAM 

tltl. 

•RTR  ROOT/*/ 

•  ISAM 

•  •99. 

•OTA  ttA**.R#T.1.0L .AtUft/IM  . 1H. . IM« . 1*1 , IR*/ 

•  ISAM 

•  •99. 

c 

•  ISAM 

•  •99. 

Mill 

•  ISAM 

1*11. 

1 •#?*•# 

•  ISAM 

•«9t. 

tAU.IT.I*)  MTI  II 

•  ISRM 

t*N. 

MRIftlROOT. «#|  • 

•  ISRM 

••  A*RM#T| !•*, *9RM#ti  tilt  SA  '.!*,'  I*  TO#  tMAtt 

■  ISRM 

•  •Mr 

•  *  )RTfla#at ATtM  1  •  AtOTI V  .  *  ) 

•  ISRM 

•  •99. 

••  ?•  1## 

HIM 

till. 

■  •  CMTIRVI 

•  1  BRM 

•  •M. 

N#l TC| ROOT, 191  COOT 

■  ISRM 

•  •If. 

IS  A9RMATI INI ,S9*. *••*,/) 

•  ISRM 

•••9. 

c 

•  ISAM 

••M. 

«  IU»(  f  tttim  IT  ■•••lit  IMCII  >  Ml  IITIMIIITI  II  IITIIK 

C  1  OAM 

•  •Mr 

•  ISAM 

Ill* . 

C 

4I44M 

1*4  1 

9«C*  I9INI  •  B|  1  )  1/49 .0 

9I44M 

IMt 

44  14  1*1.99 

414AM 

••«a 

T 1 |) «l| 11*9144911- II494C 

4I44M 

IM4  . 

14 

C9NTINUI 

4I44M 

1t*4  . 

C4UL  9T449II. T.T.TI.N.A9) 

4  14AM 

1444  . 

C 

BI44M 

t|«t. 

c 

INITIAkllf  4k  t 

9  14  AM 

till  . 

c 

4I44M 

1444  . 

4Ni |MM* 11/2 

9  ISAM 

1444  . 

44  24  U*  1  .Mas 

4  ISAM 

•  ••1  . 

24 

*.(■<>•••' 

91  BAM 

*41  Tl IN4MT . 29)  NAMI.MAMV. | 4L 1 J ) . J» 1 . MM) 

9  1  BAM 

•  ••a . 

29 

.•avail/.  l»«  .»*.»«  .»*/ )•>,  1*1  t»-l.««.iaiat> 

4  I  9  AM 

itt# . 

09  30  J «  1  ,  MM 

•  I44M 

1044  . 

20 

Akl  Jl atlANR 

4I44M 

(Ill  . 

Ai | i 1 *9i 

41  BAM 

IM1  . 

Ail  MM) *4  k 

9  1  4  AM 

1444  . 

c 

914AM 

IIM 

c 

MND  MAI  444  MIN  : 

919AM 

1440. 

c 

• 

B  19  AM 

CALL  at42lVI.A9.VM49.IN4) 

414AM 

(••a. 

CALL  MINIVI . 49. VMIN. IN4) 

914AM 

t««a . 

A  V • 1  .1*1 VM41 - VMI N ) 

919AM 

I9I4V.LT . 1 .9-24)  1 49T V • 1 

914AM 

IMS  . 

c 

914AM 

III!  . 

c 

HOT  s 

•  19AM 

144V  . 

c 

•  I  BAM 

mi . 

•  ISAM 

1444  . 

IMI  49  T  1 .10.1)  49  TO  29 

914AM 

tan . 

Cl*  1  VI  |«4)*VMIN)/4V 

419AM 

nit . 

C  1  • 2 • • 1 C 1  -  .91 

•  14AM 

aava 

00  TO  2T 

914AM 

t«ia 

34 

Cl  *o. 

919AM 

1914  . 

17 

■  ■  ON* 1 C 1  *  1  .  1*2.4 

■  19AM 

iaia . 

4k 1*1*2 

419AM 

ina . 

■  14AM 

tan . 

DO  29  1 • 1 ,K 

414AM 

1474  - 

39 

4UI  1*2 

mi . 

39 

C4NT INUC 

419AM 

mo. 

MITII.UT  .  J.  1  TIJl.TIIJl.tailll.I.I.HHI 

39 

.aava.l  ia«,.ie.*.ii,»»*,l«.iaiail 

414AM 

iui . 

4k  1 4  1  * IL4NK 

tan . 

I9| I49T .90.1)  04  TO  4# 

414AM 

1444  . 

00  41  1*2.4 

4I44M 

im . 

4  1 

4i  1  1  )  *4L4N4 

•  14AM 

1444  . 

40 

C4NTINW4 

•  14AM 

1441  . 

DO  40  1 • 1 .MM 

•  14AM 

1444  . 

90 

4LI  1  I • OOT 

•  14AM 

1444  . 

4  L (  1  ) *9LU4 

•  14AM 

1444  . 

Ait  MM ) *9104 

914AM 

W4ITIIN0UT . 40)  { AL 1 1 ) . 1 ■ 1 .MM) 

414AM 

144  2  . 

90 

.aiMii  1  iom.isi  m- ) .  »«  .  iatai ) 

•  14AM 

1441  . 

c 

•  14AM 

1444. 

c 

•  ISAM 

1444  . 

BltAM 

1444  . 

WRIT! INOUT . TO)  VMIN.VM4S 

•  ISAM 

1447  . 

70 

9  04MAT I29K.910.4. TOS .910.4) 

B  1 4  AM 

1444 

ioo 

COOT INUI 

•  14AM 

1444  . 

MKTW44 

•  19AM 

1404. 

■  40 

•  14AM 

1401  . 

tuaaouTiaa  mim.i.  i  v  .  a. aru* .  c»»«  . can  i  ,  aaam.TaaKt  .  la.TI 

9  1 AAM 

•  ISAM 

1442. 

■  ISAM 

•  ISAM 

1444  . 

4 

414AM 

1444  . 

■  ISAM 

1447  . 

INNUT:  1  *  V9CT0N  04  K  VALUI4.  L 4  NOTH  * N 

•  ISAM 

1444  . 

•  14AM 

1444  . 

IT  •  NOWS  ALkOCATOO  TO  V  IN  COikINO  44994AM 

■  ISAM 

N  •  NUM4I4  04  S-V4LUI9 

•  ISAM 

1411. 

N4UN  -  NWM4B4  04  9W4CTI449  TO  44  9L4TT90  l<*4) 

•  ISAM 

Chan  ♦  vccton  or  lcnbtn  t  contairino  cmaoactkos  to  if 

BltAM 

1411. 

U9«0  IN  T Hi  4L0T.  |CHA4|1)  19  U9C0  404  THt  4I44T 

•  ISAM 

FUNCTION.  CHAN  12)  404  TNI  BtCQNO.  BTC.) 

•  ISAM 

•  ISAM 

14*4  . 

CANT  1  ON  A90VI  THC  NkOT  12044  444MAT ) 

419AM 

KNAMC.TNAMI  •  2  VICTB49  ON  ktOBTN  20  TO  99  U9I9  494 

■  14AM 

LA49LS  ON  TNI  R  AND  V  AKI4  12041  444MAT) 

BltAM 

1  ON T  -  lONTtO  * • >  OMAkiBBT  NOINT  WILL  99  OlOBCTkV  OH  4219 

•  14AM 

1  ON  T » 1  SMALL9ST  NOINT  Wllk  94  BklSMTLV  AWAT 

•  ISAM 

414AM 

1422 

•  ISAM 

1421 . 

SU9N0UTINI9  CALLKO.  MINM4I 

•  ISAM 

1424  . 

•  ISAM 

1424  . 

N40CNAMMI4:  4N|L  SNICT04 

•  ISAM 

1424  . 

•  14AM 

1427  . 

•  ■  • 

•  14AM 

1924  . 

•  ISAM 

1924  . 

4MN4MKI30) . YN4MII20) . VN4M|1<44) .1|N) . V| IV.N4UN) 

•  ISAM 

1410. 

DATA  CANI . NLUS , 0A4H. SM. 4LANN/1HI . IN*. 1N« . 1HM, »H  / 

•  14AM 

1921  . 

00  1  1*1.19 

•  ISAM 

1912. 

1 

VNAM21 I  1  ) • 01  AUK 

•  ISAM 

1921. 

00  2  1*1.20 

•  ISAM 

1914  . 

2 

MMaMBlI  1B« 1  1 •»»*«* I  1  1 

•  ISAM 

1924  . 

DO  2  1*1.14 

•  ISAM 

1424. 

2 

•  ISAM 

1427  . 

kN*  74 

•  14AM 

1914  . 

k  V  •  44 

•  ISAM 

1419 

kVN1*kV*« 

•  14AM 

1944. 

LVM1 *LV  -  1 

•  ISAM 

1441  . 

CALL  MINMA2IV, 1 V . N . N4UN . VMIN . VMAS ) 

•  ISAM 

1442. 

can  Mianaao.a.a,  i  .  avia. aval  I 

•  14AM 

1942. 

INI  IMT.  94.  1  )  VMI  4*  VMI  4-  I  .  1*||  VMAS-VMlOI/NkOATIkV-  1  )  | 

•  19AM 

1944. 

I4| | Off. 94. 1 ) NMI 4 • NMf  4  *  1  .  1 • | | HMA 2 • NMt 4 1 / 4 k 0 AT f L H • 1  )  » 
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OUTPUT:  PaS.PSO'BTS  *  LRW8A  OUARTILB.  MBP1AN.  ANP  UPPIR  PUARTIL 
RBAR.PD  •  SAMPLB  MB AN  ANO  STANOANO  OBVIATION 

PUPPRPPRAMA  CALLBO:  OM I CK , OBBTAT . PPLPT . POPNC . WSPACI . RAO . PCOOBA 


COMMON  /PANM/  OSTAP.  OBTAIN 

01 MBNA  I  ON  NAMBISO) . I PORMt SO) . L I CAPS ( 90  ) 

0 IMS NS  I ON  S|NOIM|.VtS*t).ON|OSS|,lHS«a).POISSS).OISBai. 

P  «« t ( Alt ).««!( BAP |  , WPS (AAt) 

0  ININS  I  ON  L A0( Pol. LAO  I  ISO) .LAOS! SO) . L APP I  SO) . LAP 1P| BO ) 
DIMENSION  NAMUIB) .NAMIOia) .NAMCWSia) . ICASI(S) 

0 IMP NS I ON  10UA0T|A|.«Pt0IPI.BPU|P|.STRTiai . SNOf S t . CMANL 1 8 ) 

PATA  BPU  /  P.PI .P.PA.P. 1P.P.1A .P.TA .P.AP.P.AA.P.AA  / 

OATA  i ICAPB/ANNORM. ANAL  , AMP IPO , AMNTAL , ANL OP I . 4 MPT  1C , 

B  ANOOL  . ANBIP  . AMMO  A . ANSC I P . AMCAUC . ANMV 

0  ANSI  V . ANAL MB . ANL  PA  , AMNSAM, RNPARB , ONTO 

0  ANNS  IP  ,  A  NOLL  .  4NNAL  P  .  AN  LOS/ 

OATA  BlR/‘  ’/ 

OATA  StST.SNO.CNANL  /-. S . O ... S. O .,*•*.***  / 

OATA  AOTSN/tN«/ 

BATA  LAOt/' INPO* . ‘RMAT* . *  I  VS  • .  'OMAN * , • T I LI  * ,  •  '.!«•'  */ 

OATA  LAOO/ANO*  •,!•«  .AN  AT  ,  ANN  •  ,  S*  AM  .AN.  0. 

♦AN*  *  .SPAN  .AN  AT  .ANN  •  , V*  AN  / 

OATA  iAOIN/ANNNN  .AN  OIS . ANTNIS . ANNT 10. ANN  Bf.ANU) 

♦  0* AN  .AN  CAS.ANt  .  0»AN  / 

OATA  NAMM.OAMI#/  AM  .AN  N .  AN  IN.ANlM)  / 

OATA  NAMCtNO/AN  ,  ANOl  0  )  / 

OATA  OOCOAC/O. 1/ 

■UNIT  *  0 
ON  S  • *  1 , »t 
0  L  AO  1 1 I ♦A ) a  L  AOI  I  ) 

MNITSIRWRIT. tOOO I 

POO  PORMAT I///. TAO. • PULL  V  NON • PANAMA TN I C  ANAL VO  10 • /TOO, Of  I  IN* I  ) 
OOTOIN  OBACRIPTIVB  OTAT 1ST ICA  DON  MNARPNPSP  OATA 
PONS  •  1 . B*  tO 

01  *  A . • ATANl 1 . | 

TMOPI  •  0|  ♦  01 
SOI  •  1./0I 

NSOTOI  •  1./BONTITMOOI) 

N#«N-1 
NtOlsN 
00  IN  I  •  •  ,  N 
10  t|| t«P| 1 1 

COAL  ONICPIN.V) 

TMINPff If 

•ALL  010 TAT  I T , N, MAMS . LAO , LOUANT , O. I . PBA , P AO . RTA ) 


Ml  t»« .0*00 
00  ON  1*1,0 
Oft*«t*Mtl)< 


mi 

ae  4Mt  *  1 • v c  a » 

•  IS  AM 

aas* 

TIOI./II  ■  1  •■>•-»«»  1  1 

•  ISAM 

im 

44  94  1 ■ 1 .4 

•  ISAM 

t*ii 

441 1 1 •I4M| | | • **• 1 ott* 

•  ISAM 

tan 

•  ISAM 

nn 

a*  W«3tll»4N|l) 

•  ISAM 

aat* 

44  44  1*1,4 

•  ISAM 

a*ts . 

IF IWNtl  1  I  . 64  .  • 1  .  )  64T4  *4 

•  ISAM 

aat*. 

44  NR 3 1 1 1 . 

•  ISAM 

aatt . 

•0  DO  *4  1 • 1 . 4 

•  ISAM 

tat*. 

IH»«llN-  1  •  1  )  LE  .  1  .  1  60TQ  74 

•  ISAM 

aat* . 

*0  WK3I4* 1* t )• f . 

•  ISAM 

aa*o 

70  CONTINUE 

•  ISAM 

**•» . 

C4t  4  9*10T{U ,«K>. 44*1 . ■LK.4.  . 4.  ,1.0.414. 3, -I..1 

. 3.STRT .END. 

•  ISAM 

ta*a 

«  CHRRL  . N4MU , NAM14 . NAME .  LAB  1  ) 

•  ISAM 

aa*a 

007*  1  •  1 . • 

•  ISAM 

aas* 

4 1 1 4 1  1 )  •  4F 1 NOI ON .N4P 1 . EKUl 1 ) ) 

•  ISAM 

aat*. 

7*  eONTINUt 

•  ISAM 

aa*« . 

WRITE (NUNI T , SO)  EIU.CBI4 

•  ISAM 

aa*t. 

•  0  F0RM47I//T 10.  ‘  V  '  . 31  .  *1  FS  .  *.  1* )  , 

•  ISAM 

aa** 

•  /TIO. ' IOlU) ‘ . SB.SIFS.B. 1K| .///) 

•  ISAM 

aat* . 

•  ISAM 

aa**. 

1CA*4| * |«1 >CA*|| RASB-1  | 

•  ISAM 

aa*  i 

I  CAS* I  3 ) * l 1 CAICIKA4I ) 

•  ISAM 

aa*a . 

C 

COMPUTE  AMO  HOT  RAW  SPACING*  1*4)  A40  RAW  F©  1*1 

./•) 

•  ISAM 

aa*a 

CAll  OTOPOlON . U.NOP 1 ,WK1 .SPCFAC 1 

•  ISAM 

zat« 

c 

COMPUTE  ANO  PLOT  WEIGHTED  SPACING*  FOR  CASt  ICASC 

•  ISAM 

aa*s . 

DO  1 OO  I  •  1 , NO 

•  ISAM 

aa** . 

WR3II)  *  POPNC! Ul |« 11 . 1NUL | 

•  ISAM 

aa*t . 

1 OO  CONTINUt 

•  ISAM 

aa** 

CALL  NlPOC I 1  Nil  , 0 , NO  P 1  ,  WK  1  . WK3 . U . SW) 

•  ISAM 

aa*«  . 

c 

PLOT  CUMULATIVE  WII6HTI0  SPACING*  WITH  D ♦  AND  O' 

•  ISAM 

>400  . 

CALL  RSO 1  0 .  U .  NOP  1  .  DM. UM. OP  . UP ) 

•  ISAM 

1401  . 

CALL  FCODlAl DP . *Hj  F7 . 4 )  .144*13)) 

•  ISAM 

>402 

CALL  FCODlAl UP . SHIFT  4 )  .L*0*(*)) 

•  ISAM 

7*01 

CALL  FCODCAI  DM,  4M{F7  . 4  )  ,LA*t(  10)1 

•  ISAM 

2404 

CALL  F COOCA I UM.  *H  |  FT  .  4  )  .LA0*|I4)) 

•  ISAM 

a«o« . 

L A414I4 I • ICASC ( 1 ) 

•  ISAM 

7*0*  . 

1*11*11  )  *  1 C  AS  C ( 2 ) 

•  ISAM 

2401  . 

CALL  OPLOTIU.O. NOP 1.4LK.0. .0..1.0.4LK.3.0.1..I 

,  © .  .  1  .  . ASTER . 

•  ISAM 

a*os . 

4  NAMU . MAMCMS  .  L AS IO.  LAI*  ) 

•  ISAM 

**©• . 

NIITIIHUNIT . 1 30) 

•  ISAM 

3410. 

130  FORMAT!////// | 

•  ISA(4 

a*  1 1 

•  ISAM 

a«ia. 

CND 

•  ISAM 

a*  ia . 

S  U0R4UT INI  OUICRIH.T) 

•  ISAM 

a* «« . 

•  ISAM 

a*  n 

OU  ICR  SORT.  THIS  ALCORITH  IS  ALSO  RCFCRRCD  T|  AS  A 

PARTITIONED 

•  ISAM 

341*. 

CRCHAMG*  SORT  .  CXPCCTCO  RUNTIME  IS  PROPORTIONAL  TO 

N«L0G3!N) 

•  ISAM 

a*  1 1 . 

ALTHOUGH  TNI  WORST  CASK  IS  PROPORTIONAL  TO  ■  ••*’. 

•  ISAM 

axis. 

Rf  FERENCS  DONA  L  0  C.  RRUTM  ■  THE  ART  OP  COMPUTER  PROGRAMMING  VOL  3. 

•  ISAM 

a«  i* . 

INPUT  : 

•  ISAM 

a«ao. 

X.N  VECTOR  TO  *|  SORTED  OF  LINCTM  N 

•  ISAM 

a«a i . 

OUTPUT  : 

•  ISAM 

>422. 

•  ISAM 

a«aa. 

SU4RDUTINCS  CALLED  :  NONE 

•  ISAM 

a***. 

■  ISAM 

*43*  . 

REAL  TINI.T 

•  ISAM 

a«a*. 

INTEGER  IP,  LV|  IS)  .  1 V|  IS)  ,LP.  IUP 

•  ISAM 

a*at. 

t»! 11*1 

•  ISAM 

3421 

IVf  1  ).N 

•  ISAM 

a«a* 

IP*  1 

•  ISAM 

a«ao. 

IO  IF |  IP . LT  ,  1  )  GO  TO  7* 

■  ISAM 

a«a  t . 

IS  IF(  ( IVl  IP  I -LV|  IP )  )  . LT .  1  I  60  TO  30 

■  ISAM 

3*33  . 

GO  TO  a* 

■  ISAM 

*43* 

30  1P.1P-I 

•  ISAM 

3434  . 

CD  TO  IO 

■  ISAM 

343* 

3*  L  P  *  L  V 1  IP  )  -  i 

•  ISAM 

343*. 

1  UP  *  I V |  IP) 

•  ISAM 

a*at . 

T«T| IUP ) 

■  ISAM 

343*. 

30  IF!  |  IUP-1P )  .  LT  .  3  )  GO  TO  4* 

■  ISAM 

343*  . 

LPiiPAI 

■  ISAM 

344*. 

IMTIlPI.ll.1l  GO  TO  30 

■  ISAM 

344  4 

T  |  I  UP  )  •  T  1  L  P  I 

■  ISAM 

3443  . 

3*  IF | ( IUP-LP). LT. 3)  GO  TO  40 

■  ISAM 

3443 

1  UP • 1  UP  -  1 

■  ISAM 

3444  . 

1 F 1 T | IUP  )  . GC  . V  )  SO  TO  3* 

■  ISAM 

344*  . 

T l L  P 1 • T I  1  UP  ) 

•  ISAM 

344*  . 

GO  TO  30 

■  ISAM 

3447 

40  lUPtlUP-f 

■  ISAM 

*44* 

4*  T |  1  UP  I • T 

•  ISAM 

344* 

I  F !  (  IUP-LVI  IP  )  )  . LT .  I  I v(  IP ) - IUP )  )  GO  TO  SS 

•  ISAM 

34*0. 

GO  TO  *4 

• 

■  ISAM 

34*  1  . 

SS  LV« IPA1 ) ■ L V| IP) 

■  ISAM 

34*3  . 

IVl 1 P ♦ 1 ) • IUP* 1 

34*3. 

L  V (  IP)*  IUP*! 

■  ISAM 

34*4  . 

60  TO  70 

•  ISAM 

34**. 

40  LV| IP*1 ) • IUP*1 

■  ISAM 

34*4. 

|V( IPAI I • | V( IP) 

■  ISAM 

34*7  . 

IV! |P)| IUP* 1 

•  ISAM 

3***. 

70  |Pa|P«| 

34«*  . 

GO  TO  14 

•  ISAM 

34*0. 

7*  RETURN 

■  ISAM 

344  1  . 

404 

•  ISAM 

34*3. 

■  ISAM 

3443. 

c 

•  ISAM 

*4*4. 

■  ISAM 

*44*. 

c 

•  ISAM 

a***. 

C 

0  IN 

•  ISAM 

344t  . 

C 

THE  N- VECTOR  RR  . 

■  ISAM 

344*  . 

c 

•  ISAM 

3*4*. 

c 

•  ISAM 

34**. 

c 

•  ISAM 

3431  . 

c 

SU4R0UTINK  CALLED  j  4RD3 

•  ISAM 

*433. 

c 

•  ISAM 

*473 

•  ISAM 

*434. 

c 

■  ISAM 

343*. 

DIMENSION  K|N)  ,IR(N| .DM  1*04,9)  .W3I444.3) . T 1 *40  I 

■  ISAM 

343*. 

c 

■  ISAM 

*437. 

c 

■  ISAM 

*43*. 

c 

■  ISAM 

*43*. 

c 

■  ISAM 

*••• . 

04  10  1 ■ 1 , N 

■  ISAM 

*4*1  . 

•  11  1. 1 »•*! 1  I 

■  ISAM 

*44*. 

•  1  1  1 .*1 ■FLOAT!  1  ) 

•  ISAM 

»*•*. 

10  C4NTINU* 

•  ISAM 

*4*4. 

cal l  ••■aiwi .*,*••) 

■  ISAM 

»4**. 

c 

•  ISAM 

*4*4. 

c 

■  ISAM 

*443. 

r 

■  ISAM 

9444. 

c 

•  ISAM 

*444. 

c 

•  ISAM 

*444. 

•  ISAM 

•  441  . 

•31  I,  *  t  RON*  |  1,  1) 

■  ISAM 

*44*. 

•til , *  1 • F  fc  OAT tl| 

■  ISAM 

*44*. 

■  ISAM 

*444. 

c 

•  ISAM 

*444. 

c 

■  T 

■  ISAM 

9444. 

c 

•  ISAM 

*443. 

c 

■  ISAM 

9444. 

1*  * 

■  ISAM 

*4*4. 

JO  R • 1 • 1 

•  ISAM 

( 


I 


■  BAL  OITNM1IC  TIltaMfil  PAC1M:  illUMI  1*0  IP  AM  I  >  I  .  1C  .  •  .  I  •  • 

WT||T*AISf«»|t .  I  I ••SIR.  II  I 

I  a | NT II T .At .!,•*•)  SI  Tl  II 
«!■«  I 

II  R»*R-I 

rs*r*i 

uma. it.  •  )  ii  n  •• 

SUM*!. 

II  II 

IVMHVIMllU.t) 

•I  CIMTIIWI 

•I  tl  J»l ,  AS 

TUidUM/riMiuil 

tl  CUT  HUB 
•I  Tl  •• 

•I  »(  1  MUM  ,  B) 

•I  1*1 

IBM  .  IT  .1)  II  Tl  SI 

wTiiTt*  ah  (wan.  1 1  •warn*  i  .in 

|P|WTBSTt  It. 1.1*11  II  TO  til 
V|N)*Pl*«*  Ml 

CAB  ATI  VK  Til  II  CINTAININ6  IANRS  CORAICTID  Ml  TIBI 
til  ••  Ml  1*1,1 

K«M»t  |  I  .  S  1*0.1 

IKMliTll  ) 

1 «•  CIMTIIWI 

BIO 

IUIAIU1IIC  IPIMWtN.IMl.lUMO) 


INPUT:  IAMRI  •  TlB  VICTOR  IP  I AMR  f  OP  Til  INPBPBWDIWT  VAIIAIlB 
I  •  Tit  IWMIIR  IP  PA 1 18 0  OISCIVATIONI 
•  WTPWT:  INI  •  IPB AIMAI *  0  IAMR  COIIBIATtOM  CIBPPIC1BMT. 


COMMON  /OAT Al/  1 1 III )  ,  V  t III I  , I AIR 1 1 III ) . IANRT I  III ) 


SUMS* I . 

DO  II  1 ■ 1 .1 

DIP*  IAMR  HI  I  MIAMI  V(  I  ) 

IUMO  *  I  IIMO*0  IP*  1 1 P 
II  CIWTIMUC 

COMPUTB  SPI AIMAI *  I  IAIR  CORRELATION  COIPPICIINT 
*IWM| 

SDIMOM.P  LOATH).  |PtOAT|M-N  M  I  .  I 
INO*  t  . •SHtfM/IIBMIM 
II  TUN 
BID 

SUOMlUT 1 MB  TI|M|B,V. RMSO. VMCO . ROBL . N . NI*N ) 


•Wl PI III AM  TO  TRIM  A  llTARIATf  BATA  BIT  IP  AT  MOST 
ROBL  •BVTIBMI*  POUTS  BASSO  OR  DISTANCI  PROM  TIP 
MSB  1  AN  IN  TNI  ■  All  V  OUBCTIONS  OUT. 

INPUT:  S.V  -  BATA  IP  SIS!  N 

RMSO.TMBI  -  MID  IANS  OP  I  All  T 

ROIL  *  MAR  I MUM  MUNICH  OP  POINTS  TO  OILSTB  PI 
OUTPUT:  R.T  -  TIIMMBO  BATA  BP  Sill  NlWN 
BUI PI 06 I AMI  CALLII:  MAM . Ml N . 00 1 CR . OIOS 


BIMINI I  ON  MIN) .V(N| .IIIM(I«,|| .MBSt (SI) . 10BL I  HI . W (till 
NBMN • N 

ip  I RBBt  . i i  . •)  ■■ Twin 

RCNR • I 

1 P I Kill . HI  .  J  )  II  Tl  s 
RCNR ■ I 
ROB  L ■ a 

S  IPIRBBL.LT. N)  IB  TO  I 

MR  I Tl 1 1,1) 

I  P BIMAT | TIM . ‘IBB L  IB  IIBATII  THAN  OR  BIUAL  Tl  N'./, 

*111.  'ROIL  NAB  OBCN  BBT  BIUAL  Tl  A.  'I 
RBBL •« 

I  IP(RBIL.eT.ll|  ROBL • SO 
RMOO*MOO(ROCL .1) 

IPIRMOB.ll.fi  ROIL • RBI L  •  1 

R0JM«RBIL 

KHALP»ROBL/a 

BB  II  1*1, RNALP 

J*N- 1*1 

IBLMI  I  .  I  )ilMIB*im 
BB  L Ml  I . J ) • PLlAT ( I ) 

OlLMl ROBL *  1*1 ,  I  )  *■( Jl-RMBO 
BBLMIRBBL • I ♦ » . t ) • P t BAT I J ) 

II  CINTINBI 

CALL  BIBS | OB LM, ROBL , II) 

BB  SB  1*1, RNALP 

IBBLI 1 )• 1PIRI IBLMI RNALP* I ,!)*•.■) 

>•  CONTINUI 

I P I RCNR . B| , l|  II  Tl  l| 

II  SI  1*1.1 
VV(lt'VtS) 

SI  CUT  I  NUB 

CALL  MIBIV.N. VMIN, IMIN) 

CALL  M4B|*.I. VMAB , I MAI) 

•I  Al  1*1. RNALP 
IBLMI  I . I  )  ■ TMBB  *  TMI I 
BBLMl  1.9) • PLlAT I IMIN) 

V | IMINI • VMAB 

CALL  MINlV.N, VMIN. IMIN) 

A#  CUT  HUB 

II  II  I  •  I ,  N 
Tlimwil) 

II  CUT  HUB 

CALL  MINI  V  ,N, VMIN, IMIN) 

U  10  t* I, RNALP 

IILM(  ROBL  *|Alal|* VMAB • VMSB 

IBiMf  RML  •  I*  I  ,  1  )  tPLBATI  IMAB) 

VI  |MAI)*TMII 

CALL  MAS I  V , N, VMAB , IMAB ) 

•I  CANT  f  NUB 

II  II  I • t . I 


ISAM 

HAM 

SIAM 

HAM 

HAM 


ISAM 

HAM 

HAM 

HAM 

HAM 

HAM 

HAM 


HAM 


HAM 

HAM 

HAM 


ISAM 

HAM 


HAM 

HAM 

HAM 

HAM 

HAM 

HAM 

HAM 

HAM 

ISAM 

ISAM 

HAM 

HAM 

HAM 


HAM 
HAM 
HAM 
HAM 
HAM 
HAM 
HAM 
HAM 
HAM 
HAM 
HAM 
HAM 
HAM 
HAM 
HAM 
HAM 
HAM 
HAM 
HAM 
HAM 
HAM 
HAM 
HAM 
HAM 
HAM 
HAM 
HAM 
HAM 
HAM 
HAM 
HAM 
HAM 
HAM 
HAM 
HAM 
HAM 
HAM 
HAM 
I  SAM 
HAM 
HAM 
HAM 


HAM 

HAM 

HAM 

HAM 

HAM 

HAM 

HAM 

HAM 


ISAM 

HAM 

HAM 

HAM 

HAM 

HAM 

HAM 

HAM 

HAM 

HAM 

HAM 

HAM 

ISAM 

HAM 

HAM 

HAM 

HAM 

HAM 

HAM 

HAM 

HAM 

HAM 

HAM 

HAM 

HAM 

HAM 

HAM 

HAM 

HAM 

HAM 

HAM 

HAM 

HAM 

HAM 

HAM 

HAM 

HAM 

HAM 


o  vm»wti) 

CALL  •MI((llM.K||i,IA| 

I  OIL  I KNAL F •  I  I 

o  comti nui 

RN1lRNftir«l 
DO  ••  III, KMALP 
DO  M  J*KN1.K0IM 

IP|  (DDL  I  I  I  .Ml  .  IDfi  |  Jl  |  DO  TD  DO 

IDUUhNO 

KOI l • ROIL  *  > 

O  COMTINUI 

00  01  1*1  .  KOI** 

3  ROIL  I  I } • FLOAT tlOli(ll) 

CALL  OU1CK(KOIM.  MOIL  I 
DO  OA  1*1  .KOI** 

14  IDIU  I  |  •  IM  I  I  IDti  (11*0  .  «  I 
It  COMTINUI 

IMKCNK.IO.il  KOI  t  •  1 
MM  1  •  N  •  1 

00  too  1*1, KOli 
INO* IDIL IKOIL • l« 1  ) 

NIWNiNIWN- 1 
IMIND.IO.Nl  CO  TO  03 
00  ••  J • IMD.MM1 
R|J|*K|J*1) 

)0  COMTINUI 
13  I|MIWN*1 1*0.0 
V  I  MINN ♦ 11*0.0 
MM  1  ■  MIRIM 
>0  COMTINUI 

WRIT* l A.  I  10) 

0  FORMAT |  IOX .  *TN*  FOLLOWING  FOIMTt  WIRI  OILITIO  FROM  TMI  ORTA  OIT s  *  I 
WRITIIA.130)  I  I  OIL  I  I),  1*1, K OIL) 

IO  FORMAT  I  10K , All  I 

WRITIIA.130)  KDIL.NIWM 

IO  FORMAT t/,10R,ll,*  FOIMTS  WIRI  OILITIO  LCAVlNO  '.10,'  FOIMTt  *. 

♦  * IN  TNI  DATA  OIT. * ) 

RETURN 
■  NO 

OUOROUT I  Ml  WOFACI I Wit . CWKt , MOF 1 . MO , FOOO . U . 0  I IA, AVLWX ) 


1U0  ROUT  INI  TO  COMFUTI  OIU).  CUMULATIVI  0*0.  AMO  SIOMAO 
FOR  TN|  MOOIL  01 U) «MU»I 1 OMA-OO(U) 

IMFUT  : 

MO. MO  ?  VICTOR  OF  LlMOTM  HO  CONTAINING  LITTLO  OlU) 
FOOO  :  MVROTNItltlO  OONOITV  OUANTILI  FUNCTION 
U  *  VICTOR  OF  LlMOTM  NO  COMTAININO  U  VALUII 

OUTPUT  : 

Wit  :  VICTOR  OF  LlMOTM  MO  CONTAINING  O(U) 

CWMO  :  VICTOR  OF  LlNGTM  NO  CONTAINING  TNI 
CUMULATIVI  0*0 

0100  :  COMPUTIO  VALUI  OF  OIOMAO  •  CWMtlNO) 

OUOROUT I  Nit  CAL  L 10  :  MONl 


01 MR MS  I  ON  FOOOIMOP1 )  . MttMOPI I .UIM0P1 ) . WM 0  I  MOP  1 t , CWMO I  NOP  1 ) 
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